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ABSTRACT 


This  report  documents  the  work  conducted  in  phase  I of  an  ongoing  research 
program  of  which  the  objective  is  to  develop  a head-neck  simulator  with 
omni-directional  biofidelity.  The  report  is  composed  of  two  volumes: 

Vol.  I : Analysis  volunteer  tests 
Vol.  II:  Mathematical  simulations 

In  Vol.  I a detailed  analysis  is  presented  of  a large  number  of  human  volun- 
teer tests  conducted  by  the  Naval  Biodynamics  Laboratory  in  New  Orleans. 
These  human  subjects  were  exposed  to  frontal,  lateral  and  oblique  impacts 
with  an  impact  severity  up  to  15  g and  17  m/s.  This  analysis  results  for 
each  impact  direction  in  a simple  analog  system  that  completely  specifies 
the  observed  dynamical  behaviour. 

Vol.  II  documents  the  validation  of  the  proposed  analog  systems  and  compares 
the  behaviour  of  the  Part  572  and  Hybrid  III  head-neck  systems  with  the 
human  volunteer  behaviour.  It  follows  that  both  neck  designs  are  stiffer 
than  any  of  the  volunteers  in  the  impact  ranges  tested.  Preliminary  results 
are  presented  of  computer  aided  design  activities  in  order  to  improve 
existing  designs  in  view  of  these  new  findings. 
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CHAPTER  1 


INTRODUCTION 


1 . 1 Project  description  and  objectives 

The  National  Highway  Traffic  Safety  Administration  (NHTSA)  has  large  quan- 
tities of  test  data  on  the  head  and  neck  response  of  human  volunteer  sub- 
jects in  frontal,  oblique  and  lateral  sled  tests.  The  Vehicle  Research  and 
Test  Center  (VRTC)  will  use  these  data  for  the  development  of  an  improved 
head-neck  motion  simulator.  Where  possible,  these  human  test  data  will  be 
augmented  by  data  of  human  cadavers  tested  at  exposure  levels  beyond  those 
generally  acceptable  for  human  volunteers. 

One  or  more  candidate  mechanical  systems  will  be  identified  to  reproduce 
observed  head-neck  motions.  Currently  identified  models  will  form  the  basis 
for  this  analysis  but  the  study  will  not  necessarily  be  limited  to  such 
models.  The  performance  of  candidate  mechanical  neck  simulators  will  be 
analyzed  and  optimized  using  mathematical  model  simulations  resulting 
finally  in  design  specification  for  an  improved  dummy  neck. 

This  report  presents  the  results  of  phase  I of  this  ongoing  research  pro- 
gram. The  report  is  composed  of  two  volumes: 

Volume  I : Analysis  volunteer  tests 
Volume  II  : Mathematical  simulations. 

In  Vol.  I an  analysis  is  presented  of  a number  of  human  volunteer  tests  in 
three  different  impact  directions  conducted  at  the  Naval  Biodynamics  Labo- 
ratory (NBDL)*.  This  analysis  results  for  each  impact  direction  in  an 
analog  system  that  completely  specifies  the  observed  dynamical  behaviour. 
Vol.  II  documents  the  validation  of  the  proposed  analog  systems  by  means  of 
mathematical  simulations  and  presents  a number  of  preliminary  computer  aided 
design  activities  carried  out  in  order  to  improve  existing  dummy  necks  in 
view  of  these  new  findings. 


* NBDL:  Naval  Biodynamics  Laboratory  in  New  Orleans,  previously  Naval  Aero- 
space Medical  Research  Laboratory  Detachment  (NAMRLD) 
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1 . 2 Siinwnary  research  program:  mathematical  simulatioas 

During  this  study  numerous  mathematical  simulations  were  conducted.  The  most 
interesting  ones  will  be  summarized  here.  These  calculations  were  conducted 
with  the  MADYMO  Crash  Victim  Simulation  program  package  (1). 

Chapter  2 presents  simulations  performed  for  validation  of  the  analog  sys- 
tems proposed  in  Vol.  I.  Three  2-pivot  linkage  systems  have  been  proposed 
in  Vol.  I:  one  for  the  behaviour  in  frontal  impacts,  one  for  oblique  im- 
pacts and  one  for  lateral  ones.  Simulations  will  be  presented  here  for  all 
three  impact  directions.  Results  of  some  of  the  parameter  variations  con- 
ducted in  order  to  improve  the  insight  in  the  system's  behaviour  will  also 
be  presented  in  this  chapter. 

Chapter  3 deals  with  a model  for  the  Part  572  dummy  head-neck  assembly  in 
frontal  impacts.  The  neck  section  of  this  model  consists  of  7 segments. 
Model  results  will  be  compared  with  the  Part  572  pendulum  calibration  test 
and  with  sled  tests. 

In  chapter  4 the  behaviour  of  the  Part  572  and  Hybrid  III  dummy  head-neck 
assemblies  will  be  compared  with  the  human  volunteer  behaviour.  Significant 
differences  between  dummy  and  volunteer  can  be  observed.  The  model  of  the 
Part  572  dummy  presented  in  chapter  3 will  be  used  to  analyze  the  effect  of 
possible  design  changes. 

The  major  findings  of  these  model  simulations  will  be  summarized  and  dis- 
cussed in  chapter  5.  Based  on  these  findings  recommendations  are  made  for 
improvement  of  the  procedure  used  to  analyze  human  volunteer  test  data. 
Also  a preliminary  design  recommendation  will  be  made  for  a dummy  neck  with 
frontal  biofidelity. 
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CHAPTER  2 


SIMULATION  OF  THE  HEAD-NECK  MOTION  OF  HUMAN 
VOLUNTEERS  BY  A TWO-PIVOT  MECHANISM 


2 . 1 Introduction 

In  this  chapter  the  validation  of  the  system  for  the  human  head-neck  beha- 
viour proposed  in  Vol.  I will  be  presented.  The  analysis  is  limited  to  the 
loading  phase  of  the  motion.  A summary  of  the  system  parameters  of  the  two- 
pivot  mechanism  specified  in  Vol.  I will  be  presented  in  section  2.2,  while 
the  representation  of  this  mechanism  in  the  MADYMO  Crash  Victim  Simulator 
is  described  in  section  2.3.  In  section  2.4  results  of  the  validation 
efforts  are  presented.  An  analysis  of  deviations  between  model  and 
experiments  is  given  in  section  2.5.  Using  the  model,  the  effect  of  a number 
of  parameter  variations  was  evaluated.  Some  of  the  most  interesting  varia- 
tions will  be  presented  in  section  2.6. 

2 . 2 Summary  of  the  system  parameters  provided  in  Vol.  I 

Table  1 summarizes  the  most  important  linkage  parameters  of  the  analog 
systems  proposed  in  Vol.  I (Fig.  1).  The  orientation  of  the  head  principal 
axes  and  the  location  of  occipital  condyles  and  center  of  gravity  will  be 
specified  in  agreement  with  the  assumptions  made  in  Vol.  I: 

- the  angle  between  principal  z-axis  and  head  anatomical  z-axis  is  -36° 
(principal  z-axis  rotated  backward) 

- occipital  condyles  relative  to  head  anatomical  coordinate  system: 

x=-l.lcm  y=0cm  z=-2.6cm 

- center  of  gravity  relative  to  head  anatomical  coordinate  system: 


x = 1.2  cm 


y = 0 cm 


z = 2.9  cm 
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LmkoQe  mechanism 


thrust  vector  of 
accelerator 

Fig.  1 Two-pivot  model  for  the  human  head-neck  system. 

Fig.  2 illustrates  the  location  of  the  head  center  of  gravity  relative  to  a 
coordinate  system  oriented  along  the  principal  axes  with  the  origin  near  the 
occipital  condyles. 

The  mass  and  moments  of  inertia  of  the  head  in  the  model  will  be  based  on 
estimated  values  for  subjects  H00083  and  H00093  (see  Table  6,  Vol.  I).  The 
mean  value  of  these  data  (including  instrumentation)  is  used  as  model  input 
and  is  incorporated  in  Table  1. 
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center  of  rotation 
(LOWER  JOINT  I 


frontal 


Fig.  2 


Orientation  of  coordinate  systems 
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Table  1 Linkage  parameters  for  the  loading  phase  based  on  30  impact  tests 
with  two  subjects  (H00083  and  H00093). 


Parameter 

Frontal 

Lateral 

Oblique 

direction 

direction 

direction 

Geometrical  parameters 


upper  pivot  location 

occipital 

occipital 

occipital 

condyles 

condyles 

condyles 

link  length 

0.125  m 

0.125  m 

0.125  m 

lower  pivot  location 

see  table  8 

see  table  8 

see  table  8 

(Vol.  I) 

(Vol.  I) 

(Vol.  I) 

initial  head  link  flexion 

0 

0 

0 

initial  head  link  torsion 

0 

0 

0 

initial  neck  link  rotation 
Joint  parameters 

14  deg. 

5 deg. 

10  deg. 

lower  pivot  linear 

stiffness 

upper  pivot  linear 

1.25  Nm/deg. 

1.5  Nm/deg. 

1.25  Nm/deg. 

torsion  stiffness 
upper  pivot  free  range 
of  motion  in  the  plane 

0.25  Nm/deg. 

0.5  Nm/deg. 

of  impact* 

Mass  distribution 

head  center  of  gravity 

mass 

I 

jXX 

lYY 

zz 

25° 

10° 

see  fig.  2 
4.28  kg 
0.0230  kg. m2 
0.0238  kg. m2 
0.0124  kg. m2 

o 

o 

* Joint  is  locked  if  maximum  excursion  is  reached. 


2 . 3 MADYMO  model  set-up 

The  model  consists  of  three  elements  (torso,  neck  and  head)  connected  by  two 
ball  and  socket  joints.  The  element  and  absolute  coordinate  systems  are 
illustrated  in  Figure  2.  For  the  lateral  and  oblique  simulations  the  torso, 
neck  and  head  links  are  rotated  90  and  45  degrees,  respectively,  about  the 
z-axis  of  the  absolute  (laboratory)  reference  frame.  The  initial  orientation 
of  the  segments  in  the  present  MADYMO  version  (i.e.  version  3)  has  to  be 
specified  by  so-called  Bryant  angles.  For  the  lateral  and  oblique  simula- 
tions these  angles  were  obtained  by  a pre-processing  program. 
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For  the  simulations  presented  here  a new  joint  model  was  developed:  the 

flexion-torsion  joint  model  (Fig.  3).  Consider  the  joint  between  element  i 

and  element  j.  In  both  elements  vectors  t and  f are  defined.  By  means  of 

these  vectors  it  is  possible  to  specify  two  rotations  Of  and  3,  respectively. 

The  rotation  a is  the  angle  between  the  vectors  t.  and  t.  and  is  called  the 

J 

flexion  angle.  This  angle  is  defined  in  the  interval  from  0 to  180  degrees. 

As  a consequence  of  the  rotation  a the  vector  f.  takes  a new  position  f.'. 

-J 

The  angle  p is  defined  as  the  rotation  around  the  t.-axis,  which  transforms 

f.'  into  f..  This  angle  is  called  the  torsion  angle  and  is  defined  from  -180 
-J  -J 

to  + 180  degrees.  In  this  way  it  is  possible  to  specify  two  torques,  i.e.  a 
flexion  and  a torsion  torque. 


Fig.  4 shows  the  t and  f vectors  used  for  the  human  head-neck  system.  A more 
detailed  description  of  this  joint  model  is  presented  in  Appendix  A. 
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Fig.  4 Definitions  of  the  t and  f vectors  in  the  two-pivot  model. 

This  new  joint  model  can  be  used  for  the  lower  pivot  using  the  lower  pivot 
linear  flexion  stiffness  specified  in  Table  1 as  torque-flexion  characteris- 
tic. Torsion  in  the  lower  pivot  will  be  suppressed  since  the  upper  pivot 
already  allows  a torsion  motion. 

For  the  upper  pivot  the  situation  is  more  complicated  due  to  the  required 

locking  behaviour.  For  the  loading  phase  of  this  joint,  i.e.  the  phase  where 

the  angle  a increases,  in  a standard  way  the  torque-rotation  function  as 

illustrated  in  Fig.  5a  can  be  prescribed.  The  free  range  of  motion  is 

different  for  each  impact  direction  as  illustrated  in  Table  1.  The  joint 

stop  stiffness  is  arbitrarily  selected  as  500  Nm/rad.  For  unloading  (i.e.  a 

becomes  smaller)  the  same  characteristic  can  be  used  until  a is  reached.  If 

o 

a becomes  smaller  than  the  original  torque-rotation  relation  has  to  be 

adjusted  in  order  to  obtain  the  locking  effect.  This  new  characteristic  is 

illustrated  in  Fig.  5b.  In  some  of  the  mathematical  simulations,  which  will 

be  treated  later,  it  was  found  that  the  joint  stop  of  the  upper  pivot  was 

not  reached  and  consequently  the  locking  mechanism  was  not  activated.  Due  to 

this  the  joint  model  was  modified  in  such  a way  that  the  locking  mechanism 

is  activated  even  if  a doesn't  reach  to  Of  -value.  Torsion  characteristics 

o 

for  the  upper  pivot  are  defined  according  to  Table  1. 
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flexion 

torque 


Fig.  5 Flexion  characteristic  of  the  upper  pivot. 

For  the  flexion  motion  of  the  upper  pivot  damping  will  be  specified  which  is 
activated  in  the  joint  stop.  To  investigate  the  influence  of  this  damping 
two  simulations  of  a frontal  volunteer  test  (representative  for  LX3544, 
LX3548  and  LX3550)  were  conducted: 

- DAMP  1 : undamped  model 

- DAMP  2:  model  with  a damping  coefficient  of  1 Nms/rad. 

Significant  oscillations  are  present  in  the  accelerations  in  case  of  absence 
of  damping  (Fig.  6).  The  effect  of  this  damping  on  the  rotations  of  the  neck 
and  head  element,  however,  is  small. 


Fig.  6 


Model  results  from  a simulation  of  a frontal  volunteer  test  for 
an  undamped  (DAMPl)  and  a damped  upper  pivot  joint  (DAMP2). 


In  the  flexion-torsion  joint  model  a third  direction  of  motion  can  be  dis- 
tinguished, namely  perpendicular  to  the  flexion  direction  (Fig.  7).  Mathe- 
matical simulations  showed,  that  after  the  locking  mechanism  is  activated, 
the  tj  vectors  tended  to  slide  away  in  the  third  motion  direction  along  a 
cone  (Fig.  7).  As  a consequence,  out  of  plane  motion  will  occur  and  the  t^ 
vector  can  move  to  the  other  side  of  the  cone  dependent  on  the  loads  acting 
on  the  head-neck  system.  To  prevent  this  undesirable  motion,  linear  damping 
was  specified  in  this  third  direction.  For  the  simulations  described  in  this 
chapter  a damping  coefficient  of  10  Nms/rad  is  selected. 


Fig.  7 Third  motion  direction  in  the  flexion-model. 


2 . 4 Validation  of  the  two  pivot  model 


The  model  runs  for  the  three  impact  directions  with  data-sets  as  described 
in  the  preceding  sections  will  be  called  the  reference  simulations.  In 
Appendix  B the  input  data-sets  are  presented.  These  simulations  will  be 
compared  with  a set  of  tests  under  similar  conditions: 


- frontal  impacts:  LX3544,  LX3548  and  LX3550 

- lateral  impacts:  LX1831,  LX2013  and  LX2355 

- oblique  impacts:  LX3133,  LX3153,  LX3122  and  LX3158. 

All  these  tests  are  of  a moderate  severity.  The  test  conditions  are  presen- 
ted in  Table  3 in  Vol.  I.  As  part  of  the  sensitivity  study  (section  2.6), 
the  most  severe  volunteer  test  (test  LX3616)  in  the  present  data  base  will 
be  simulated  separately.  An  approximation  for  the  T1  acceleration  measured 
in  the  volunteer  tests  will  be  used  as  model  input  (Fig.  8). 


The  following  output  parameters  will  be  used  for  model  validation: 

- head  torsion  and  neck  link  rotation  as  function  of  the  head  flexion, 

- head  and  neck  rotations  as  functions  of  time  as  defined  in  Fig.  1, 

- linear  accelerations  of  the  head  anatomical  origin  with  respect  to  the 
laboratory  reference  frame  as  function  of  time, 

- head  angular  accelerations  as  function  of  time. 


Model  results  (REFl)  and  test  results  for  the  frontal  impact  situation  are 
presented  in  Figs.  9-10.  Fig.  9 shows  that  head  and  neck  rotation  functions 
are  predicted  very  well  by  the  model.  A small  difference  can  be  observed  for 
the  head  flexion  0 and  the  corresponding  angular  accelerations,  which  start 
too  early  in  the  model.  Due  to  this  the  maximum  relative 
angle  between  the  head  and  neck  element  is  only  16  degrees,  while  the 
volunteer  tests  showed  an  average  value  of  25  degrees. 


degrees , 
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T1  accelerations  and  model  input  in  frontal,  lateral  and  oblique 
direction. 


Fig.  8 
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TIME  IN  MSEC 

Fig.  9 Rotations  for  the  frontal  volunteer  tests  and  the  reference  model 
run  (REF1=DAMP2) . 


u 


Fig.  10  Linear  accelerations  of  the  anatomical  origin  in  laboratory  direc- 
tions and  head  angular  accelerations  for  the  frontal  volunteer 
tests  and  the  reference  model  run  (REF1=DAMP2) . 
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Results  of  the  simulation  (REF2)  for  the  lateral  volunteer  tests  appear  to 
be  slightly  closer  to  the  volunteer  tests  than  the  results  of  the  frontal 
test  simulation  (Figs.  11  and  12).  The  same  holds  for  the  simulation  of  the 
oblique  tests  (REF3)  with  one  exception:  the  head  torsion  tjj  is  predicted 
by  the  model  is  too  small  (Figs.  13  and  14).  The  maximum  relative  flexion 
angle  in  the  upper  pivot  appears  to  be  predicted  adequately  by  both  the 
lateral  and  oblique  model  set-ups  (Figs.  11  and  13). 
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Fig.  11  Rotations  for  the  lateral  volunteer  tests  and  the  reference  model 
run  (REF2). 
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Fig.  12  Linear  accelerations  for  the  anatomical  origin  in  laboratory 
directions  and  head  angular  accelerations  of  the  lateral  volunteer 
tests  and  the  reference  model  run  (REF2) . 
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Fig.  13  Rotations  for  the  oblique  volunteer  tests  and  the  reference  model 
run  (REF3). 
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Fig.  14  Linear  accelerations  of  the  anatomical  origin  in  laboratory  direc- 
tions and  head  angular  accelerations  for  the  oblique  volunteer 
tests  and  the  reference  model  run  (REF3). 
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2 . 5 Conclusioas  from  the  validation  study 

The  preceding  simulations  show  two  significant  deficiencies  in  the  model: 

- the  initial  maximum  relative  upper  pivot  rotation  in  the  impact  plane  is 
too  small  in  case  of  frontal  impacts, 

- the  head  torsion  p is  not  predicted  well  by  the  model  in  case  of  oblique 
impacts . 

For  the  first  deficiency  the  following  explanation  can  be  given:  the  initial 
joint  characteristic  of  the  lower  pivot  is  expected  to  have  a dominant  in- 
fluence on  the  initial  relative  rotation  in  the  upper  pivot.  The  joint 
characteristics  used  as  model  input  for  the  lower  pivot  are  based  on  esti- 
mations from  the  loads  calculated  with  respect  to  T1  rather  than  the  torso 
center  of  rotation  (chapter  6 of  Vol.  I).  As  a consequence  the  selected 
joint  characteristic  for  the  lower  pivot  might  not  be  completely  correct. 
The  influence  of  this  stiffness  will  be  illustrated  in  section  2.6.2. 

The  second  deficiency  is  probably  caused  by  the  limited  accuracy  of  the 
graphical  technique  used  for  the  estimation  of  the  upper  pivot  torsion 
stiffness.  The  effect  of  a different  upper  torsion  stiffness  will  be 
illustrated  in  section  2.6.4.  The  different  ways  in  which  the  head  torsion 
is  specified  in  model  and  volunteer  tests  could  also  contribute  to  the 
observed  differences.  For  the  volunteer  tests,  head  torsion  is  derived  from 
the  third  Euler  angle.  In  the  model,  the  torsion  angle  is  determined 
according  to  the  equations  presented  in  section  2.3.  The  possible  influence 
of  this  is  not  clear  yet  and  is  subject  to  further  study. 

2 . 6 Sensitivity  study  with  the  two-pivot  model 

A number  of  parameter  variations  were  conducted  in  order  to  get  a better 
understanding  of  the  system  behaviour  and  to  explain  differences  between 
model  and  experiment.  Some  of  the  most  interesting  variations  will  be  con- 
sidered here: 

- the  model  behaviour  for  a single  severe  test  using  subject  specific  data 

(2.6.1), 

- the  effect  of  lower  pivot  stiffness  (section  2.6.2), 
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- the  effect  of  the  absence  of  a locking  mechanism  in  the  upper  pivot 
(section  2.6.3) , 

- the  effect  of  the  torsion  stiffness  in  the  upper  pivot  (section  2.6.4). 

2.6.1  Simulation  of  a single  severe  test  using  subject  specific  data 

Simulations  were  conducted  to  evaluate  the  model  validity  for  a higher  im- 
pact severity  and  to  analyse  the  effect  of  subject  specific  input  data.  Two 
simulations  will  be  presented  here: 

- SENl:  simulation  of  the  most  severe  volunteer  test  (subject  H00093  in 

test  LX3616)  using  subject  specific  data  (see  Table  2), 

- SEN2:  simulation  of  the  most  severe  volunteer  test  using  the  reference 

data  set. 

In  both  simulations  the  initial  head  and  neck  link  position  was  based  on  the 
actual  position  in  the  human  volunteer  test  (see  Table  2). 

Table  2 Differences  in  system  parameters  relative  to  the  reference  simula- 
tion in  simulations  SENl  and  SEN2 . 


Parameter 

SENl 

SEN2 

- initial  head  link  flexion  in  plane  of 
impact 

-12  [deg] 

-12  [deg] 

- initial  neck  link  rotation  in  plane  of 
impact 

20  [deg] 

20  [deg] 

- lower  pivot  linear  flexion  stiffness 

1.88  [Nm/deg] 

- head  mass 

- moments  of  inertia  of  the  head: 

4.03  [kg] 

I 

XX 

0.0209 

I 

yy 

0.0215  [kg. m2] 

I 

zz 

0.0110 

Fig.  15  shows  the  T1  acceleration  of  test  LX3616  together  with  the  model 
input.  Model  and  test  results  are  presented  in  Figs.  16  and  17.  Both  simula- 
tions appear  to  predict  the  results  of  this  more  severe  test  situation  quite 
realistically.  Fig.  16  shows  that  the  model  with  the  subject  specific  data 
performs  slightly  better  than  the  model  with  the  average  linkage  parameters. 
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Fig.  15  T1  acceleration  for  the  most  severe  volunteer  test  (LX3616)  and 
model  input. 
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Fig.  16  Rotations  for  the  most  severe  frontal  volunteer  test  and  model 
runs  (SENl  and  SEN2). 
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Fig.  17  Linear  accelerations  of  the  anatomical  origin  in  laboratory  direc- 
tions and  head  angular  accelerations  for  the  most  severe  frontal 
volunteer  test  and  model  runs  (SENl  and  SEN2) . 
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2.6.2  The  effect  of  the  flexion  characteristic  of  the  lower  pivot 

Since  the  lower  pivot  stiffness  was  estimated  rather  than  calculated,  it  was 
decided  to  perform  a sensitivity  study  for  the  initial  part  of  the  torque- 
rotation  characteristic  (for  larger  rotation  values  the  torque  estimations 
are  expected  to  be  more  accurate).  One  of  the  changes  (SENS)  is  shown  in 
Fig.  18. 


- 9 - 00  /■  rad  I 


Fig.  18  Flexion  characteristic  of  the  lower  joint  in  the  reference  model 
and  in  model  SENS. 


From  a comparison  between  Figs.  9 and  19  it  can  be  concluded,  that  the  model 
results  are  improved  with  respect  to  the  maximum  relative  angle  in  the  upper 
pivot.  In  these  simulations,  the  relative  angle  reaches  2S  degrees  rather 
than  the  16  degrees  achieved  in  the  reference  simulation.  Based  on  this 
relatively  large  influence  it  is  recommended  to  calculate  lower  neck  torques 
directly  at  the  torso  center  of  rotation  rather  than  at  Tl.  Other  model 
results  were  hardly  effected  by  this  variation. 
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Fig.  19  Rotations  for  the  frontal  volunteer  tests  and  the  model  run  with 
a changed  flexion  characteristic  of  the  lower  joint  (SENS). 


2.6.3  The  effect  of  the  locking  mechanism 


One  of  the  most  interesting  findings  from  the  analysis  of  the  volunteer  test 
was  the  existence  of  a so-called  locking  mechanism  in  the  upper  neck  region. 
This  mechanism,  most  likely  due  to  muscle  activity,  was  represented  in  the 
model  by  a relatively  stiff  torque-rotation  characteristic  of  500  Nm/ rad 
(joint  stop),  which  is  activated  as  soon  as  the  upper  pivot  reaches  the 
locking  or  maximum  position.  In  tests  with  human  cadavers  (no  muscle 
activity)  or  for  much  higher  impact  levels  this  locking  mechanism  might 
be  partly  or  completely  absent.  The  possible  effect  of  this  is  simulated  by 
specifying  a much  lower  locking  stiffness  (SEN4:  50  Nra/rad  rather  than  500 
Nm/rad)  and  is  illustrated  in  Fig.  20.  From  this  figure  it  can  be  concluded 
that  the  neck  rotation  0 slightly  decreases  whereas  the  head  rotation  ({) 
strongly  increases  compared  to  the  reference  simulation  (Fig.  9).  In  other 
words,  if  the  locking  mechanism  is  absent  the  head  will  tip  over. 


PHi  I HE  I A 
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Fig.  20  Rotations  for  the  frontal  volunteer  tests 
locking  mechanism  in  the  upper  point  (SEN4) . 


and  a model  without  a 
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2.6.4  The  effect  of  the  torsion  stiffness  of  the  upper  pivot  (oblique  test) 

To  illustrate  the  influence  of  the  torsion  stiffness  of  the  upper  pivot, 
results  of  a model  run  with  a less  stiff  characteristic  (SENS:  0.35  Nm/deg 
instead  of  0.5  Nm/deg)  will  be  presented.  Fig.  21  shows  that  the  model 
prediction  for  the  head  torsion  is  improved  by  this  change.  Since  the 
magnitude  of  this  variation  is  in  the  same  order  of  magnitude  as  the  ac- 
curacy of  the  graphical  technique  used  to  determine  this  parameter,  it  is 
recommended  to  specify  this  parameter  more  accurately  in  future.  Another 
aspect  is,  that  the  torsion  stiffness  has  been  defined  from  the  volunteer 
data  on  the  basis  of  a torsion  angle  specified  as  an  Euler  rotation  about 
the  head  local  z-axis.  In  the  mathematical  model  a different  definition  is 
used  for  the  torsion  angle.  The  effect  of  this  will  be  studied  in  a next 
phase  of  this  project. 


► 'P-'Po  Irad] 

0.30  0.60  0.90  1.20  ISO 


Fig.  21  Head  torsion  as  function  of  head  flexion  resulting  from  model 
simulations  and  volunteer  tests. 
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CHAPTER  3 


SIMULATION  OF 
A STANDARD 


THE  PART  572  HEAD-NECK  ASSEMBLY  IN 
CALIBRATION  TEST  AND  SLED  TESTS 


3 . 1 Introduction 

In  this  chapter  calibration  and  sled  tests  conducted  at  VRTC  (2)  with  the 
Part  572  head-neck  assembly  will  be  simulated  mathematically.  The  two-dimen- 
sional version  of  MADYMO  will  be  used  since  only  frontal  test  results  are 
available.  This  model  will  be  used  as  a design  tool  for  Computer  Aided 
Design  (CAD)  activities  (see  chapter  4). 

The  compound  pendulum  impact  apparatus  specified  in  FMVSS  208  Part  572  was 
used  for  the  calibration  tests.  A detailed  sketch  of  this  set-up  is  shown 
in  Fig.  22.  Instrumentation  includes  a three-accelerometer  block  located 
at  the  center  of  gravity  of  the  head  and  a fourth  accelerometer  mounted  on 
the  pendulum  itself  (see  Fig.  22  for  exact  location).  Pendulum  velocity 
shortly  before  impact  is  measured  by  a light  trap  mounted  at  the  accelero- 
meter location.  Head  motion  is  measured  with  two  rotary  potentiometers,  one 
mounted  at  the  pendulum  base  and  the  other  located  at  the  head  c.g..  The 
rotary  potentiometers  are  connected  by  a linear  potentiometer  (Fig.  23). 

The  sled  tests  were  conducted  on  a HYGE  sled  using  a specially  designed 
buck  for  mounting  the  head-neck  assembly  (Fig.  24).  The  impact  severity 
for  the  sled  test  is  much  higher  than  the  severity  of  the  calibration  test, 
22  m/s,  35  g compared  to  11  m/s,  22  g.  The  same  instrumentation  for 
measuring  head  motion  and  accelerations  was  used  as  in  the  neck  pendulum 
calibration  tests  (Fig.  23).  The  tests  are  described  more  in  detail  in  (2). 

In  the  next  section  a description  of  the  model  set-up  is  presented.  A more 
detailed  description  can  be  found  in  (3) . Model  verification  will  be 
presented  in  section  3.3. 
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Fig.  22  Test  set-up  of  the  standard  calibration  test  for  the  Part  572 
head-neck  assembly. 


Fig.  23  Potentiometers  for 
measuring  the  head-neck  motions. 
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3.2  Model  description 

The  basic  set-up  of  the  model  is  shown  in  Fig.  25.  Element  1 represents  the 
sled  or  pendulum.  Elements  2 through  7 represent  the  neck,  whereas  the  head 
is  represented  by  element  8.  Model  input  is  the  displacement  of  the  neck 
base  (element  1),  i.e.  the  displacement  (acceleration)  of  pendulum  and  the 
sled,  respectively.  The  pendulum  acceleration  and  velocity  were  measured  at 
0.20  m above  the  neck  base  (Fig.  22).  Model  input  acceleration  follows  from: 


a . 

L 


a 

m 


(1) 


where : 

a . = model  input  acceleration 
1 

r^  = distance  between  neck  base  and  pendulum  pivot  (1.841  m) 

r,  = distance  between  accelerometer  and  pendulum  pivot  (1.657  m) 
b 

a^  = measured  pendulum  acceleration 


General  model  set-up 


Model  specifications 
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Joint  stiffness 
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Fig.  25  General  model  set-up  and  specifications  for  the  model  of  the  Part 
572  head-neck  assembly. 
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The  pendulum  accelerations  in  4 tests  together  with  the  selected  model  input 
are  presented  in  Fig.  26.  For  the  sled  tests  a simple  approximation  of  the 
sled  acceleration  can  be  used  as  input.  Fig.  27  shows  an  envelope  of  the 
sled  accelerations  in  10  similar  sled  tests  and  the  corresponding  model 
input . 


Fig.  26  Pendulum  accelerations  and  model  input. 


Time  (msec) 


Fig.  27  Sled  accelerations  and  model  input. 

The  neck  which  is  considered  to  be  a uniform,  symmetrical  rubber  cylinder  is 
divided  into  6 segments  of  equal  length  (Fig.  25;  E2-E7).  The  total  length 
of  the  neck  is  0.111  m (excluding  the  steel  endplates).  The  total  neck  mass 
is  divided  and  assigned  to  the  6 rigid  neck  elements  and  partly  to  the  head 
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and  torso  element.  Mass  distribution  data  for  the  head  were  derived  from  a 
study  of  the  mass  distribution  of  the  complete  Part  572  dummy.  The  results 
of  the  study  are  incorporated  in  Appendix  C.  The  most  significant  model 
parameters  are  summarized  in  Table  3 and  Fig.  25. 

Table  3 Masses  and  moments  of  inertia  for  the  2D  model  of  the  Part  572 
head-neck  assembly. 


mass 

moment 

of  inertia 

[kg] 

[kg. m2] 

head 

element 

4.48 

0.0307 

neck 

element 

0.125 

0.0010 

MADYMO  allows  specification  of  the  initial  angular  velocity  of  an  element. 
For  the  Part  572  head-neck  assembly  in  the  calibration  tests  this  angular 
velocity  can  be  obtained  from: 


(2) 


where : 

a»Q  = the  initial  angular  velocity 

Vq  = the  initial  linear  velocity  at  the  accelerometer  location  (7.42  m/s) 
r^  = the  distance  between  pendulum  pivot  and  the  accelerometer  (1.841  m) 

The  resulting  angular  velocity  (i.e.  4.03  rad/s)  is  applied  to  the  elements 
representing  the  head  and  neck  (elements  2-8).  In  addition  gravity  is 
prescribed  in  the  model. 

The  joint  model  in  MADYMO  2D  allows  specification  of  elastic  joint  torques, 
viscous  damping  and  Coulomb  friction.  An  expression  for  the  elastic  joint 
torque-rotation  characteristic  in  each  joint  of  the  model  can  be  obtained  by 
treating  each  segment  of  the  neck  as  an  uniform  linear  elastic  bar  loaded 
with  moments  and  concentrated  forces  at  each  end  (Fig.  28).  If  small  dis- 
placements are  assumed  for  each  element  the  deflection  equation  for  element 
i yields: 
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Ay.  = 


SI 


(Ax)3 


(Ax)' 


3EI 


2EI 


(3) 


A(|).  = 
1 


SI 


(Ax)^  . Ax 


2EI 


El 


Y 

I 


X 


Fig.  28  Bar  loaded  with  concentrated  forces  and  moments. 


(4) 


Parameters  in  the  above  relations  are  explained  in  Fig.  28.  The  flexibility 
of  this  elastic  segment  is  represented  in  the  model  by  a joint  located  in 
the  segment  center  (Fig.  29).  An  elastic  torque  function  is  defined  for  this 
joint  as  follows: 

M = k . (j) 

Then  the  rotation  of  the  joint  can  be  described  by: 


A(l)*  _ ^i  ^si 
i ~ 2k. 


. Ax 


(5) 
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Fig.  29  Model  elastic  bar. 
Similarly,  from  Fig.  29: 

* Ax  * 

Ay  . = — sin  (A(J) . ) 

1 2 1 


For  small  displacements  sin 


(A0.)  = A0. 
1 1 


so , 


1 


(6) 


(7) 


To  determine  a value  of  which  makes  the  linkage  behave  like  an  elastic 
bar,  equate  the  expressions  obtained  for  rotations  ((4)  and  (5)).  This 
yields : 


* 


SI 


(Ax)' 


“bi 


Ax 


2«bi 


+ F 


SI 


Ax 


2EI 


El 


2k. 

1 


This  reduces  to: 


k. 

1 


Ax 


(8) 
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Alternately,  equate  the  expressions  for  displacements  ((3))  and  (7)).  This 
yields : 

AY.  = AY? 

1 1 


F . . (Ax)^  . . (Ax)^  F . . (Ax)^  M,  . . (Ax) 

SI  D1  _ SI  D1 

3EI  2EI  " 4k.  2k. 

1 1 


or: 


k. 

1 


Ax 


3F 


SI 


4F 


SI 


Ax 


Ax  + 


If  Ax  is  quite  small  this  relation  reduces  to: 


k. 

1 


Ax 


(9) 


This  is  identical  to  the  expression  obtained  from  consideration  of  rota- 
tions. An  average  value  of  EI=8.32  [Nm^]  was  determined  from  published 
static  bending  test  data  (4,5).  The  corresponding  elastic  joint  stiffness  is 
523  [Nm/rad] . The  damping  coefficient  required  for  critical  damping  in  each 
joint  was  estimated  by  analyzing  a one-degree-of  freedom  system.  This 
analysis  resulted  in  a critical  damping  coefficient  for  flexion  of  1.45 
[Nms/rad] . 

3.3  Model  validation 


In  order  to  validate  this  model,  corridors  have  been  determined  using 
results  from  the  sled  tests  and  the  calibration  tests.  The  following  para- 
meters will  be  used  for  validation: 

- the  head  rotation, 

- the  horizontal  (x)  and  vertical  (z)  displacements  of  the  head  c.g.  rela- 
tive to  the  neck  base, 

- the  horizontal  (x)  and  vertical  (z)  accelerations  of  the  head  c.g. 
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Three  different  simulations  have  been  conducted  for  the  calibration  tests, 
namely  with  undamped,  "50%  critically"  (0.75  [Nms/rad])  and  "critically" 
damped  (1.45  [Nms/rad])  neck  joints.  Results  from  these  simulations  are 
shown  in  Fig.  30.  Reliable  model  results  are  obtained  for  the  head  rota- 
tions, the  vertical  head  displacement  and  the  linear  accelerations.  Rela- 
tively large  deviations  can  be  observed,  however,  for  the  horizontal 
displacements,  which  can  be  explained  by  the  absence  of  shear  and  elonga- 
tion in  the  model.  From  the  model  results  it  appears  that  the  "50%  cri- 
tically" damped  model  provides  the  most  realistic  results.  So  this  model 
will  be  used  for  further  validation  at  a higher  impact  severity  level  (sled 
test).  Fig.  31  shows  the  results  of  this  simulation  together  with  the 
experimental  corridors.  Here,  too,  quite  realistic  results  are  obtained 
from  the  simulations,  except,  for  the  horizontal  displacements,  as  could 
be  expected. 
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head  rotation 


corridors  of  test  results 

undamped  model  (Part  1) 

— "50%  critically"  damped  model 
(Part  2) 

— "critically"  damped  model 
(Part  3) 


X-  displacement 


z-displacement 


Fig.  30  Model  and  test  results  of  the  standard  calibration  tests. 
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Fig.  31  Model  and  test  results  of  the  sled  tests 
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Based  on  these  simulations  it  can  be  concluded  that  this  model  offers  a 
realistic  simulation  for  the  behaviour  of  the  Part  572  dummy  head-neck 
assembly  in  frontal  impacts.  Although  the  model  significantly  underestimates 
the  horizontal  head  displacement,  this  is  not  considered  to  be  a major 
drawback  for  use  of  the  model  as  a computer  aided  design  tool. 
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CHAPTER  4 


COMPUTER  AIDED  DESIGN  OF  AN  IMPROVED  MECHANICAL  HEAD-NECK  ASSEMBLY 


4 . 1 Introduction 

The  Part  572  (Hybrid  II)  neck  was  designed  to  perform  realistically  in 
frontal  impacts.  In  a study  by  Muzzy  and  Lustick  (6)  displacements  of  the 
Hybrid  II  bead-neck  system  were  compared  with  human  volunteers  in  10  g 
frontal  impacts.  The  complete  dummy  was  tested  in  a restraint  system  simi- 
lar to  that  of  the  human  volunteers.  Significant  differences  in  the  dis- 
placements of  the  head  center  of  gravity  and  the  head  angle  relative  to  the 
torso  were  observed:  the  Hybrid  II  neck  was  found  to  be  stiff er  than  the 
volunteer ' s . 

In  section  4.2  a more  direct  comparison  between  the  Part  572  head-neck 
assembly  and  human  volunteer  behaviour  will  be  presented  by  prescribing 
the  measured  human  volunteer  T1  acceleration  directly  to  the  base  of  the 
neck.  In  this  way  test  data  are  not  diffused,  for  instance  by  the  thorax 
compliance.  The  more  recent  Hybrid  III  head-neck  system  which,  like  the 
Part  572  head-neck  assembly,  was  designed  for  frontal  biofidelity  is  also 
incorporated  in  this  experimental  evaluation. 

Section  4.3  deals  with  preliminary  computer  aided  design  activities  con- 
ducted in  this  study.  Results  will  be  presented  of  design  changes  intro- 
duced in  the  Part  572  neck  in  order  to  improve  the  mechanical  behaviour  in 
view  of  the  findings  from  the  volunteer  tests. 

4 . 2 Comparison  between  existing  neck  designs  and  volunteer  necks 

The  Part  572  and  Hybrid  III  neck  were  both  designed  to  perform  realistic- 
ally in  frontal  impacts.  At  the  Vehicle  Research  and  Test  Center  an  expe- 
rimental study  has  been  conducted  as  part  of  the  "New  Car  Assessment  Test 
Procedure  Analysis"  project,  in  which  both  necks  were  exposed  to  various 
impact  levels  (2).  In  one  of  the  tests  the  sled  acceleration  pulse  was 
based  on  the  most  severe  NBDL  human  volunteer  test,  i.e.  subject  H00093  in 
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test  LX3616.  Both  dummy  head-neck  assemblies  were  equipped  with  one  trans- 
lational and  two  rotary  potentiometers  in  order  to  determine  the  relative 
head  displacements. 

Fig.  32  presents  occipital  condyle  trajectories  calculated  from  the  instru- 
mentation data  together  with  the  human  volunteer  trajectory  of  test  LX3616. 


Fig.  32  Occipital  condyle  trajectories  and  head  rotation  in  two  dummy  neck 
designs  and  the  most  severe  human  volunteer  test. 


For  the  dummy  head-neck  assemblies  the  occipital  condyles  are  assumed  to 
coincide  with  the  center  of  the  head-neck  connection.  Head  rotation-time 
histories  for  dummies  and  volunteer  are  also  shown  in  this  figure.  It 
follows  that  head  rotation  and  horizontal  o.c.  displacements  are  slightly 
smaller  for  both  head-neck  assemblies  than  they  are  for  volunteers.  Much 
larger  deviations,  however,  can  be  observed  for  the  vertical  displacements: 
both  neck  designs  appear  to  be  stiff  in  comparison  with  the  volunteer 
behaviour*.  These  differences  between  dummy  and  volunteer  are  even  more 
apparent  if  one  considers  that  subject  H00093  shows  relatively  small  head 
excursions  compared  to  subject  H00083. 


* In  ref.  (7)  and  ref.  (8)  center  of  gravity  trajectories  for  dummy  and  vo- 
lunteer were  compared  rather  than  the  occipital  condyle  trajectories. 
Due  to  an  error  in  the  data  processing,  however,  volunteer  c.g.  traject- 
ories and  peak  c.g.  displacements  were  not  presented  completely  correct 
in  these  references.  The  effect  of  this  in  view  of  the  large  deviations 
between  dummy  and  volunteer  behaviour  is  small. 
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4.3  Design  variations 


By  means  of  the  mathematical  Part  572  head-neck  model  (PART2)  presented  in 
the  preceding  chapter,  the  effect  of  various  changes  in  the  Part  572  neck 
design  will  be  evaluated.  Results  of  4 interesting  design  principles  will 
be  presented  here: 

CADI:  reduced  neck  stiffness  (50%) 

CAD2:  reduced  neck  stiffness  (50%)  and  a free  range  of  motion  with  locking 

mechanism  in  the  head-neck  connection 

CAD3:  longer  neck  with  a reduced  neck  stiffness  (50%)  where  the  neck  is 

connected  to  the  head  4 cm  above  the  center  of  gravity 
CAD4:  very  soft  neck  (stiffness  about  10%  of  original  Part  572  neck  stiff- 

ness) reinforced  at  the  front  and  back  by  two  relatively  stiff 
springs . 


Fig.  33  illustrates  these  design  principles  together  with  the  reference 
model.  In  order  to  introduce  the  locking  mechanism  in  the  second  variation, 
the  existing  joint  model  in  MADYMO  2D  was  modified  in  accordance  with  the 
simulation  of  the  locking  mechanism  in  the  2-pivot  model  presented  in 
section  2.3. 


Reference  model  (section  32) 
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CAD  2 
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CAD  3 
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resistive  joint 
stiffness  260Nm/rad 

t T ? 

! '‘1  i 
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resistive  joint 
stiffness  60  him /red 

r~T  resistive  joint 

I ? s tiff  ness  260  Nm/rad 

V 

Fig.  33  Design  changes  in  the  standard  Part  572  head-neck  assembly. 
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For  each  design  change  two  simulations  are  conducted,  i.e.  a simulation  of 
the  most  severe  frontal  volunteer  test  (test  LX3616)  and  a simulation  of  a 
set  of  moderate  tests  (tests  LX3544,  LX3548  and  LX3550).  Approximations  for 
the  T1  accelerations  of  the  volunteer  tests  will  be  used  as  model  input 
(Figs . 8 and  15) . 

4.4  Model  results 

Results  of  the  model  simulations  together  with  human  volunteer  test  results 
are  presented  in  Figs.  34-37.  The  upper  sections  of  these  figures  show 
the  occipital  condyle  trajectories.  For  the  human  volunteer  tests  these 
trajectories  are  presented  relative  to  the  corrected  T1  coordinate  system. 
For  the  model  simulations  these  results  represent  the  trajectories  of  the 
center  of  the  neck  top  relative  to  a coordinate  system  with  its  origin  in 
the  center  of  the  neck  base  and  with  the  same  orientation  as  the  laboratory 
coordinate  system.  The  lower  section  shows  results  for  the  head  rotations 
as  function  of  time. 
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Fig.  34  Occipital  condyle  trajectories  and  rotation-time  histories  (CADI). 
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Fig.  35  Occipital  condyle  trajectories  and  rotation-time  histories  (CAD2). 
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Fig.  36  Occipital  condyle  trajectories  and  rotation-time  histories  (CADS). 
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Fig.  37  Occipital  condyle  trajectories  and  rotation-time  histories  (CAD4) . 


V. 
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Fig.  34  shows  the  results  for  a halved  neck  stiffness  (CADI).  Although  a 
much  better  approximation  of  the  occipital  condyle  trajectories  is  obtained 
by  this  design  change,  still  the  head  excursion  appears  to  be  too  small.  On 
the  other  hand  head  rotations  are  much  larger  now  than  in  the  human  volun- 
teer tests. 

By  means  of  simulation  CAD2  an  attempt  is  made  to  reduce  the  head  rotation 
by  introduction  of  a free  range  of  motion  with  a locking  mechanism  in  the 
top  of  the  neck  as  could  be  observed  in  the  human  volunteer  tests.  With  this 
design  change  the  o.c.  trajectory  is  not  affected  while  the  head  rotation  is 
indeed  slightly  reduced  compared  to  CADI  (Fig.  35). 

In  the  preceding  design  modifications  the  geometrical  configuration  of  the 
neck  was  not  changed.  Simulation  CAD3  considers  the  effect  of  a much  longer 
neck  (almost  twice  as  long)  which  is  connected  in  the  head  above  the  center 
of  gravity.  In  addition  the  neck  stiffness  is  halved  as  in  the  variations 
CADI  and  CAD2.  This  design  change  results  in  o.c.  trajectories  similar  to 
those  resulting  from  CADI  and  CAD2  (the  o.c.  location  relative  to  the  head 
reference  frame  is  identical  to  CADI  and  CAD2).  The  head  rotation  as  func- 
tion of  time  appears  to  be  improved  in  the  initial  phase  of  the  motion  (up 
to  125  ms)  compared  to  CADI  and  CAD2.  For  larger  time  values,  however,  the 
rotation  strongly  increases  resulting  in  a much  larger  peak  head  rotation. 
Fig.  38  roughly  illustrates  this  behaviour.  Up  to  125-150  ms  the  head  c.g. 
is  located  below  the  connection  of  the  head  to  the  neck  resulting  in  a 
negative  torque  applied  to  the  head  due  to  the  impact  acceleration  forces. 
Consequently  positive  rotation  of  the  head  will  be  limited.  For  larger  time 
values,  however,  the  head  c.g.  passes  the  head-neck  connection  and  the  im- 
pact acceleration  field  which  is  still  active  will  now  cause  a positive 
torque  on  the  head. 
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acceleration 

field 


f = 10  ms 


t - 150  ms 

Fig.  38  Kinematics  for  simulation  CAD3. 

Results  of  the  final  design  change  to  be  considered  here  (i.e.  CAD4)  are 
presented  in  Fig.  36.  The  geometrical  properties  of  this  design  are  iden- 
tical to  the  Part  572  neck.  The  neck  stiffness  is  reduced  by  a factor  of 
almost  10.  At  a distance  of  2 cm  from  the  neck  centerline  two  relatively 
stiff  springs  are  located  between  the  torso  and  the  head,  one  spring  at  the 
front  side  and  one  at  the  back  side.  In  this  way  the  head-neck  system  is 
expected  to  behave  more  or  less  like  a four  bar  linkage  mechanism  where  the 
torso,  head  and  the  individual  springs  represent  one  bar.  Results  of  this 
simulation  appear  to  be  very  realistic  (Fig.  37).  Resulting  linear  and 
angular  accelerations  are  presented  together  with  volunteer  test  results 
in  Fig.  39.  These  results  appear  to  be  quite  realistic. 
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Fig.  39  Comparison  of  accelerations  of  the  anatomical  origin  in  laboratory 
directions  in  moderate  frontal  volunteer  tests  and  in  the  most 
se^re  frontal  volunteer  test  and  simulation  CAD4. 
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CHAPTER  5 


DISCUSSION  AND  CONCLUSIONS 


The  study  presented  here  summarizes  a number  of  mathematical  simulation 
activities  conducted  in  phase  I of  the  VRTC  head-neck  research  program.  The 
activities  can  be  subdivided  into  three  categories,  namely,  the  validation 
of  the  analog  system  concept  proposed  in  Vol.  I of  this  final  report,  the 
formulation  of  a mathematical  model  for  the  Part  572  dummy  head-neck  assem- 
bly and  some  preliminary  computer  aided  design  activities. 

The  model  simulations  presented  in  chapter  2 for  verification  of  the  analog 
systems  show  that  the  proposed  2-pivot  linkage  mechanisms  provide  a realis- 
tic representation  of  the  human  volunteer  behaviour.  Model  results  appear  to 
be  close  to  the  human  volunteer  test  results  with  two  exceptions:  the 
relative  rotation  in  the  upper  pivot  in  frontal  impacts  and  the  head  torsion 
in  oblique  impacts  are  found  to  be  too  small. 

The  sensitivity  study  presented  in  section  2.6  reveals  that  an  improvement 
in  the  simulations  with  respect  to  the  above  deficiencies  might  be  expected 
if  the  lower  pivot  torque  data  (i.e.  in  the  torso  center  of  rotation)  are 
determined  directly  from  the  human  volunteer  data  rather  than  estimated  from 
T1  origin  torque  calculations.  Instead  of  linear  approximations  for  the 
pivot  stiffness  characteristics  the  use  of  piece-wise  linear  or  higher  order 
functions  should  be  considered.  The  limited  accuracy  of  the  graphical 
techniques  used  to  approximate  the  human  volunteer  torque-rotation  charac- 
teristics could  be  improved  by  using  numerical  techniques. 

The  underestimation  by  the  model  of  the  head  torsion  in  case  of  oblique  im- 
pacts might  also  be  caused  by  the  different  ways  in  which  the  torsion  angle 
is  defined  in  the  volunteer  tests  and  the  model.  In  the  volunteer  tests 
torsion  is  defined  by  third  Euler  angle,  while  in  the  model  as  a consequnce 
of  the  introduction  of  the  flexion-torsion  model,  torsion  is  defined  in  a 
slightly  different  way  (see  section  2.3).  The  influence  of  this  is  not  com- 
pletely clear  yet  and  is  subject  of  further  research.  The  system  parameters 


for  the  analog  systems  relate  to  average  estimations  from  tests  with  two 
different  subjects.  A calculation  showed  that  model  predictions  can  be  im- 
proved slightly  by  the  introduction  of  subject  specific  data  rather  than 
average  data. 

The  second  model  activity  in  this  study  deals  with  a model  of  the  Part  572 
dummy  head-neck  assembly.  The  neck  section  in  this  model  is  represented  by  6 
rigid  segments  connected  by  pin  joints.  Model  verification  was  conducted  for 
two  different  impact  severities:  the  standard  neck  calibration  test  and  a 
more  severe  sled  test.  Realistic  model  predictions  were  obtained,  except  for 
the  horizontal  head  displacements.  Differences  in  displacements  can  be 
explained  by  the  absence  of  elongation  and  shear  in  this  model.  This 
deficiency  is  not  considered  to  be  a serious  limitation  for  use  of  the  model 
as  a computer  aided  design  tool.  If  necessary  model  predictions  can  be 
roughly  corrected  for  the  underestimation  of  this  output  parameter.  If  such 
a correction  is  introduced  in  the  simulation  results  of  the  design  examples 
in  chapter  4,  model  predictions  in  general  would  improve  with  respect  to  the 
volunteer  behaviour.  In  future  simulations,  representation  of  the  Part  572 
neck  could  be  improved  by  the  development  of  a special  flexible  element  in 
MADYMO. 

Experimental  evaluation  of  the  Part  572  and  Hybrid  III  neck  system  with 
respect  to  the  volunteer  behaviour  showed  that  both  necks  are  too  stiff, 
particularly  for  the  vertical  displacement  of  the  head  (chapter  4). 
Mathematical  simulations  have  been  conducted  (see  chapter  4)  in  order  to 
optimize  the  design  of  the  Part  572  neck  in  view  of  the  findings  from  the 
volunteer  tests.  The  first  results  of  these  computer  aided  design  efforts 
are  quite  promising.  The  best  results  up  to  now  are  obtained  from  a design 
that  consists  of  a soft  neck  reinforced  by  a relatively  stiff  structure  in 
the  front  and  the  back  section.  The  overall  size  of  this  neck  is  the  same  as 
in  the  present  Part  572  neck.  A possible  realization  of  such  a design  prin- 
ciple is  a soft  rubber  cylinder  imbedded  by  strong  fibres  or  provided  with 
one  or  more  separate  stiff  rubber  segments.  Building  one  or  more  prototypes 
according  to  this  design  principle  and  experimentally  evaluating  them  under 
test  conditions  representative  of  the  most  severe  human  volunteer  test  (i.e. 
run  LX3616,  15  g,  17  m/s)  is  strongly  recommended. 
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The  primary  aim  of  the  computer  aided  design  activities  described  here  was 
the  evaluation  of  the  feasilibity  of  the  model  as  a computer  aided  design 
tool.  Results  presented  here  relate  to  only  a limited  set  of  possible  design 
principles.  A more  systematic  approach  could  lead  to  other  attractive  design 
concepts.  Such  an  analysis  should  not  be  limited  to  the  neck  design  itself 
but  should  also  address  the  neck  location  relative  to  thorax  and  head, 
the  head  mass  distribution  and  the  head-chest  interaction.  In  the  present 
analysis,  for  instance,  human  volunteer  occipital  condyles  are  assumed 
to  coincide  with  the  center  of  the  neck  top.  In  this  way,  however,  volunteer 
head  center  of  gravity  and  dummy  head  center  of  gravity  will  not  be  in 
the  same  position  (deviation  0.032  m) . 

Results  presented  here  are  valid  only  for  relatively  low  g-levels.  One 
subject  in  the  present  database  was  exposed  to  the  most  severe  impact 
condition  i.e.  17  m/ s and  15  g.  New  volunteer  test  results  which  recently 
have  been  obtained  from  NBDL,  include  several  tests  for  this  impact  seve- 
rity level.  In  addition  cadaver  test  results  should  be  obtained  for  higher 
impact  levels.  The  locking  mechanism  observed  in  the  volunteer  tests  will 
most  likely  be  absent  in  cadaver  tests.  A simulation  which  was  conducted 
in  section  2.6.3  showed  a larger  head  rotation  in  this  situation  and  the 
so-called  tipping  over  of  the  head  link  relative  to  the  neck  link. 

It  is  expected  that  computer  aided  design  activities  in  close  interaction 
with  the  development  and  testing  of  promising  neck  prototypes  will  lead 
to  a neck  with  frontal  biofidelity.  Once  such  a neck  is  obtained  the 
analysis  can  be  extended  to  a neck  with  omni-directional  biofidelity. 
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APPENDIX  A 


FLEXION  - TORSION  MODEL 


A1  INTRODUCTION 


For  the  simulations  presented  in  chapter  2 a new  joint  model  is  used.  In 
this  appendix  a brief  description  of  this  model  is  presented. 

A2  THEORY 


Consider  the  joint  between  element  i (lower  numbered)  and  element  j (higher 
numbered)  (Fig.  Al).  In  both  elements  vectors  t en  f are  defined  with  the 
following  properties: 


1 


(Al) 


Fig.  Al  Flexion-torsion  model. 

Two  angles  a and  p will  be  defined  to 
angle  a is  defined  as  the  angle  between 
flexion  angle: 


fa-  li  “ ~j 

I 

I 


calculate  the  elastic  torques.  The 
the  vectors  t.  and  t.  and  is  called 
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a = arccos 


(A.t. ,A.t .) 
1-1  ’ j-j 


(A2) 


where  0 ^ a < 180° 

A.  . = the  (3*3)  rotation  matrices  of  element  i and  i 

The  rotation  a is  a rotation  about  the  u-axis,  whereas  the  u-vector  is 
defined  as: 


u 


A.t.  X A.t . 

ni 

A.t.  X A . t . 


(A3) 


As  a consequence  of  the  rotation  a,  the  vector  f^  takes  the  position  f ^ ' as 
indicated  in  Fig.  Al.  The  vector  f ^ ' can  be  calculated  by  the  following 
equation: 


f.'  = (A.f.,u).u  + (A . f . -u(A . f . ,u) ) . cos  a + (u  X A.f.).  sin  a 

j 1-1  - - 1-1  - 1-1  - - 1—1 


(A4) 


The  angle  p will  be  defined  as  the  rotation  around  the  t^-axis,  which  trans- 
forms f j ' into  f^  and  is  called  the  torsion  angle.  This  angle  can  be  cal- 
culated from: 


p = arccos  (f.',  A.f.)  (A5) 

-J  J~J 

The  angle  P can  be  either  negative  or  positive.  A direction  of  the  rotation 
can  be  specified,  namely, 


if  ((fj’ 


X Af.),  A.t.)  < 0 
-J  J-J 


then  p is  negative. 


The  angle  P is  defined  in  the  interval  (-180,180). 

The  elastic  torques  are  specified  as  function  of  a and  p respectively: 


= f(a)  (A6) 

= f(P) 

where : 

= elastic  torques  for  flexion, 

= elastic  torques  for  torsion. 


The  flexion- torsion  model  calculates  also  a damping  torque  according  to 
the  following  equation: 


where : 

b = damping  coefficients  in  three  directions 

V = relative  velocity  for  the  calculation  of  the  damping  moment. 


(A7) 


The  relative  velocity  y consists  of  three  mutually  perpendicular  components 


h 


' «i ^ . . ; -i»M  . j-  ._  . «r-  , ,^-,  ,.i 


cip»l.^  65  ':■  !Mf  ;r^t 

i 


__ 

^M:_  .!i« 





i/j  ^ 

^ ^ :j: 

- ''ir 

w: 

APENDIX  B:  Reprints  of  input  data  sets 


Reprint  of 


6 


input  data  of  two  pivot  model  for  frontal  impact  (DAMP2  REFl) 


KUN 

Nsr 
ST  =■ 
LBR 
NU  = 
liEO 
GEO 
GEO 
HA  = 
MA  = 
HA  = 


3 1 


1 


3 2 

0.000 
0.000 
0.000 
l.OOOE+03 
5.000E-01 
4.280E+00 


0.000 

0.000 

0.000 


0.000 
0.000 
0.125 
l.OOOE^Ol 
5.000E-03 
2.300E-02 


0.000 

0.000 

0.051 


0.000 
0.000 
0.000 
l.OOOE+01 
5.000E-03 
2.380E-02 


0.000 

0.063 

0.031 

l.OOOE+01 

l.OOOE+01 

1.205E-02 


PCC  = 

0.0000 

0.0000 

0.0000 

see  =■  0.0000 

0.0000 

0.0000 

PCC  = 

0.0000 

0.0000 

0.0000 

see  = 0.0000 

0.0000 

0.0000 

PCC  = 

0.0000 

0.2444 

0.0000 

see  = 0.0000 

0.0000 

0.0000 

PCC  = 

0.0000 

-0.6284 

0.0000 

see  = 0.0000 

0.0000 

0.0000 

TO  = 0 

.00000 

TS  = 0 

.00100 

TE  * 0.30000 

NJTQS 

= 0 

NT  = 

2 

TT  !=■ 

1.  1. 

2.  0 

.000  1. 

0.  3.  0.00 

0.00 

0.00  0 

.00  0.00  0.00 

0.000 

0.000 

1.000  -0 

.242  0.000  0. 

970  1,000  0.000 

0.000  0.970  0.000 

TT  2 = 

1.  2. 

3.  0 

.440  3. 

2.  3.  0.00 

1.00 

0.00  10 

.00  0.00  0.00 

0.000 

0.000 

1.000  0 

.766  0.000  0. 

643  1.000  0.000 

0.000  0.643  0.000 

KFSH  > 

3 RFSH  = 5 

FSH  = 

2. 

2. 

2. 

FSM  = 

0.000 

0.000 

0.000 

FSH  =■ 

0. 

0. 

0. 

FSH  => 

1.000 

1.000 

1.000 

FSH  ’ 

-72. 

-500. 

0. 

RAHPl 

= 0.100 

RAHP2  = 

1.000 

0.242 

-0.766 


ACC  =■ 
ACC  = 
ACC  = 
KFSA  = 
FSA  = 
FSA  = 
FSA  = 
FSA  = 
FSA  =■ 
FSA  = 
FSA  = 
FSA  =■ 
FSA  = 
FSA  = 
FSA  =■ 
FSA  = 
FSA  = 
FSA  = 
FSA  = 
FSA  =• 
FSA  = 
NET  = 


1 0 0 
0 0 2 
0 0 2 
2 RFSA  = 
8.000E+00 
0.000E«-00 
O.OOOE+00 
4.000E-02 
O.OOOE+00 
8.000E-02 
-7.000E+01 
9.000E-02 
-1.150E+02 
l.lOOE-01 
-9.000E+01 
1.500E-01 
-1.000E^02 
2.600E-01 
O.OOOE+00 
3.000E-01 
O.OOOE+00 
3 


17 

4.000E+00 

O.OOOE+00 

0.000E<-0O 

4.900E-02 

O.OOOE+00 

5.000E-02 

-9.810E+00 

l.OOOE+00 

-9.810E+00 

O.OOOE+00 

O.OOOE+00 

O.OOOE+00 

O.OOOE+00 

O.OOOE+00 

O.OOOE+00 

O.OOOE+00 

O.OOOE+00 


ET 

1 = 

1. 

1. 

0.075 

0.150 

0.150 

0.000 

0.000 

-0.130 

2. 

0.000 

ET 

2= 

1. 

2. 

0.050 

0,050 

0.075 

0,000 

0.000 

0.062 

2. 

0.000 

ET 

3= 

1. 

3. 

0.075 

0.075 

0,100 

0.019 

0,000 

0,020 

2. 

0.000 

NPLT 

= 

1 

PIT 

1 = 

0. 

0. 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0. 

000  0.1 

0.000  0.000  0.000  0.000 


NBT 


i O 
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Reprint  of 


input  data  of  two  pivot  model  for  lateral  impact  (REF2) 


«UN  ; 
Nsr  = 
ST  = 
L8R  = 
NU  = 
GEO  = 
GEO  = 
GEO  = 
MA  = 
MA  = 
MA  » 


1 


1 

3 1 

3 

3 2 
0.000 
0.000 
0.000 
l.OOOE+03 
5.000E-01 
4.280E+00 


0.000 

0.000 

0.000 


0.000 
0.000 
0,125 
l.OOOE+01 
5.000E-03 
2.300E-02 


0.000 

0.000 

0.051 


0.000 
0.000 
0.000 
l.OOOE+01 
5.000E-03 
2.380E-02 


0.000 

0.063 

0.031 

l.OOOE^Ol 

l.OOOE+01 

1.235E-02 


PCC  = 0.0000 

0.0000  0.0000 

see  =■ 

0.0000 

0.0000 

0.0000 

PCC  = 0.0000 

0.0000  1.5707 

see  = 

0.0000 

0.0000 

0.0000 

PCC  = 0.0000 

0.0873  1.5707 

see  =■ 

0.0000 

0.0000 

0.0000 

PCC  = 0.6283 

0.0000  1.5707 

see  « 

0.0000 

0.0000 

0.0000 

TO  = 0.00000 

TS  = O.OOlOb 

TE  = 0. 

30000 

NJTQS  = 0 

NT  ’ 2 

TT  1=  1.  1. 

2.  0.000  1. 

0.  3. 

0.00 

0.00 

0.00  0.00  0.00  0.00 

0.000 

0.000  1.000  0 

.000  0 

.087  0.' 

996  1.000  0.000  0.000  1.000  0.000 

TT  2=  1.  2. 

3,  0.175  3. 

2.  4. 

0.00 

1.00 

0.00  10.00  0.00  0,00 

KFSM 


0,000  0.000 

4 RFSM  = 5 


0.000 

1.000  0.586  -0.087  0.806  1.000  0.000  0.000  0.809  0.000  -0.588 


FSM 

=■  2. 

2. 

2. 

2. 

FSM 

= 0.000 

0.000 

0.000 

0.000 

FSM 

= 0. 

0. 

0. 

0. 

FSM 

= 1.000 

1.000 

1.000 

1.000 

FSM 

= -86. 

-500. 

0. 

-14. 

RAMPl  = 0.100 

RAMP2 

= 1.000 

ACC 

= 1 0 

0 

ACC 

= 00 

2 

ACC 

= 00 

2 

KFSA  = 
FSA  => 
FSA  = 
FSA  * 
FSA  = 
FSA  =■ 
FSA  = 
FSA  =■ 
FSA  = 
FSA  = 
FSA  = 
FSA  = 
FSA  ’ 
FSA  = 
FSA  = 
FSA  = 
FSA  = 
FSA  => 
FSA  = 
FSA  = 
FSA  = 
FSA  = 
FSA  = 
FSA  = 
NET  = 


2 RFSA  = 
l.lOOE+01 
O.OOOE+00 
O.OOOE+00 
5.200E-02 
O.OOOE+00 
8.200E-02 
-2.500E+01 
9.400E-02 
-1.050E+02 


1. 

-1. 

1. 

-4. 

1. 

-6. 

1. 

-1. 

1. 

-5. 

2. 

0. 

3 

0. 

3 


090E-01 

OOOE+02 

200E-01 

500E+01 

460E-01 

500E+01 

710E-01 

500E+01 

960E-01 

OOOE+01 

600E-01 

OOOE+00 

OOOE-01 

OOOE+00 


23 

4.000E+00 

O.OOOE+00 

O.OOOE+00 

7.900E-02 

0.000E»-00 

8.000E-02 

-9.810E+00 

l.OOOE+OO 

-9.810E+00 

O.OOOE<-O0 

O.OOOE+00 

O.OOOE<-O0 

0.000E4-00 

O.OOOE+00 

O.OOOE+00 

O.OOOEMIO 

O.OOOE+00 

0.000E4-00 

O.OOOE+00 

O.OOOE+00 

0.000E>0O 

O.OOOE+00 

O.OOOE+00 


ET 

1 = 

1. 

1. 

0.075 

0.150 

0.150 

0.000 

0.000  -0.130 

ET 

2= 

1. 

2. 

0.050 

0.050 

0.075 

0.000 

0.000  0.062 

ET 

3= 

1. 

3. 

0.075 

0.075 

0.100 

0.019 

0.000  0.020 

NPLT 

= 

1 

PLT 

1 = 

0. 

0. 

0.000 

0.000 

0.000 

0.000 

0.000  0.000 

2. 

2. 

2. 


0.000 

0.000 

0.000 


NBT  = 0 


0.000  0.000  0.000  0.000  0.000  0.000 
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Reprint  of 


input  data  of  two  pivot  model  for  oblique  impact  (REFS) 


KUN  3 
NST  = 1 

sr  = 3 1 
LB«  - 3 

NU  = 3 2 1 

UEO  = 0.000  0.000  0.000  0.000  0.000  0.000 

OEO  = 0.000  0.000  0.000  0.000  0.000  0.063 

OEO  = 0.000  0.000  0.125  0.051  0.000  0.031 

MA  = 1.000E*03  l.OOOE+01  l.OOOE+01  l.OOOE+01 

MA  » 5.000E-01  5.000E-03  5.000E-03  l.OOOE^Ol 

MA  = 4.280E+00  2.300E-02  2.380E-02  1.235E-02 


PCC 

= 0.0000  0.0000  0.0000 

see  = 

0.0000 

PCC 

= 0.0000  0.0000  0.7854 

see  = 

0.0000 

PCC 

= 0.0000  0.1746  0.7854 

see  = 

0.0000 

PCC 

= 0.4746  -0.4286  0.8906 

see  = 

0.0000 

TO  = 

0.00000  TS  =■ 

0.00100 

o 

N 

UJ 

30000 

NJTQS  = 0 

NT  = 

2 

TT 

!=■  1.  1.  2. 

0.000  1. 

0.  3. 

0.00 

0.000  0.000  1.000  -0 

.123  0 

.123  0. 

TT 

2=  1.  2.  3. 

0.350  3. 

2..  4. 

0.00 

0.000  0.000  1.000  0 

.678  -0 

.123  0. 

KFSM 

> 4 RFSM  =■ 

5 

FSM 

- 2.  2. 

2. 

2. 

FSM 

= 0.000  0.000 

0.000  0. 

000 

FSM 

= 0.  0. 

0. 

0. 

FSM 

= 1.000  1.000 

1.000  1. 

000 

FSM 

= -72.  -500. 

0.  - 

29. 

RAMPl  =■  0.100  RAMP2 

= 1.000 

ACC 

= 10  0 

ACC 

= 002 

ACC 

= 002 

ICFSA 

1 = 2 RFSA  = 17 

FSA 

= 8.000E*-00 

4.000E+00 

FSA 

= O.OOOE+00 

O.OOOE+00 

FSA 

= U.OOOE+00 

O.OOOE^OO 

FSA 

= 5.U00E-02 

6.300E-02 

FSA 

= O.OOOE+00 

O.OOOE+00 

FSA 

= 6.40OE-02 

6.400E-02 

FSA 

= -2.5OOE*01 

-9.810E+00 

FSA 

= 8.600E-02 

l.OOOE-Ol 

FSA 

= -l.lOOE+02 

-9.810E+00 

FSA 

= 9.800E-02 

O.OOOE+00 

FSA 

= -1.200E+02 

O.OOOE+00 

FSA 

= 1.320E-01 

0.00OE*-00 

FSA 

= -6.500E+01 

O.OOOE+00 

FSA 

= 2.200E-01 

O.OOOE+00 

FSA 

= O.OOOE+00 

O.OOOE+OO 

FSA 

= 3.000E-01 

O.OOOE+OO 

FSA 

= O.OOOE+00 

O.OOOE+OO 

NET 

= 3 

ET 

1=  1.  1.  0. 

075  0.250 

0.150 

0.000 

ET 

2=  1.  2.  0. 

050  0.050 

0.075 

0.000 

ET 

3=  1.  3.  0. 

075  0.075 

0.100 

0.019 

NPLT 

= 1 

PIT 

o 

o 

o 

II 

000  0.000 

0.000 

0.000 

NBT 

= 0 

0.0000  0.0000 

0.0000  0.0000 

0.0000  0.0000 

0.0000  0.0000 


0.00  0.00  0.00  0.00  0.00 

985  1.000  0.000  0.000  0.992  0.008  0.123 

1.00  0.00  10.00  0.00  0.00 

725  1.000  0.000  0.000  0.731  0.008  -0.683 


0.000 

-0.130 

2. 

0.000 

0.000 

0.062 

2. 

0.000 

0.000 

0.020 

2. 

0.000 

0.000 

0.000 

0. 

000  0.000  0.000  0.000  0.000 

0.000 
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1 


Reprint  of  input  data  of  model 


of  Part  572  (sled) 


RUN 

NST 

a 

1 

ST 

8 1 

L8ft 

8 

NU 

8 7 6 S 4 3 2 1 

liEO 

3 

0.000  0.000  0.000  0.000 

GEO 

= 

0.000  0.000  0.000  0.008 

GEO 

= 

0.000  0.016  0.000  0.008 

GEO 

= 

0.000  0.016  0.000  0.008 

GEO 

3 

0.000  0.016  0.000  0.008 

GEO 

3 

0.000  0.016  0.000  0.008 

gEO 

=> 

0.000  0.016  0.000  0.008 

GEO 

3 

0.000  0.016  0.006  0.043 

HA  = 

l.OOOE+03  l.OOOE+01 

MA  = 

1.250E-01  l.OOOE-03 

MA  = 

1.250E-01  l.OOOE-03 

MA  = 

1.250E-01  l.OOOE-03 

MA  = 

1.250E-01  l.OOOE-03 

MA  » 

1.250E-01  l.OOOE-03 

MA  =■ 

1.250E-01  l.OOOE-03 

MA  > 

4.480E+00  3.070E-02 

PCC 

a 

0.0000  see  » 0.0000 

PCC 

3 

0.0000  see  =■  0.0000  . 

PCC 

s 

0.0000  see  » 0.0000 

PCC 

3 

-0.3500  see  * 0.0000 

PCC 

3 

-0.3500  see  * 0.0000 

PCC 

3 

-0.3500  see  =■  0.0000 

PCC 

3 

-0.3500  see  * 0.0000 

PCC 

3 

-0.3500  see  ’ 0.0000 

PCC 

3 

-0.3500  see  =■  0.0000 

PCC 

3 

-0.3500  see  = 0.0000 

TO  = 

0 

.00000  TS  =■  0.00100  TE  » 0 

15000 

MAXF 

20  MAXM  = 7 

NJTQS 

> 7 

MT 

!=■ 

1.  1.  2.  2.  0.750  0.000 

0.000 

0.000 

MT 

2=- 

1.  2.  3.  1.  0.750  0.000 

0.000 

0.000 

MT 

3 = 

1.  3.  4.  1.  0.750  0.000 

0.000 

0.000 

MT 

4= 

1.  4.  5.  1.  0.750  0.000 

0.000 

0.000 

MT 

5» 

1.  5.  6.  1.  0.750  0.000 

0.000 

0.000 

MT 

6= 

1.  6.  7.  1.  0.750  0.000 

0.000 

0.000 

MT 

U 

1.  7.  8.  1.  0.750  0.000 

0.000 

0.000 

KFSM 

« 

2 RFSM  « 7 

FSM 

3 

3.000E+00  3.000E+00 

FSM 

a 

-l.OOOE+00  -1.350E+00 

FSM 

a 

5.230E+02  5.230E+02 

FSM 

3 

O.OOOE+00  -3.500E-01 

FSM 

a 

O.OOOE+00  0.000E<-00 

FSM 

3 

l.OOOE+OO  6.500E-01 

FSM 

3 

-5.230E+02  -5.230E+02 

ACC 

3 

1 0 

ACC 

3 

0 2 

ACC 

a 

0 2 

ACC 

= 

0 2 

ACC 

3 

0 2 

ACC 

3 

0 2 

ACC 

3 

0 2 

ACC 

a 

0 2 

KFSA 

» 

2 RFSA  = 17 

FSA 

3 

8.000E*00  4.000E*-00 

FSA 

O.OOOE+00  O.OOOEK)0 

FSA 

0.000E«-00  O.OOOE+00 

FSA 

3.000E-02  3.000E-02 

FSA 

O.OOOE+00  O.OOOE+00 

FSA 

7.500E-02  3.100E-02 

FSA 

-3.350E+02  -9.810E+00 

FSA 

8.900E-02  l.OOOE+00 

FSA 

3 

-3.450E*-02  -9.810E*00 

FSA 

9.800E-02  O.OOOE+00 

FSA 

-3.300E<-02  O.OOOE+00 

FSA 

1.290E-01  O.OOOE+00 

FSA 

-1.650E+02  0.000EK)0 

FSA 

a 

1.470E-01  O.OOOE+00 

FSA 

O.OOOE+00  O.OOOE+00 

FSA 

2.000E-01  O.0OOE*0O 

FSA 

O.OOOE^OO  O.OOOE+00 

NET 

0 

NBT 

0 
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Reprint  of  input  data  of  first  CAD  model  (CADI) 


■ION 

1 

NST 

a 

1 

ST  > 

1 

8 1 

L3R 

3 

3 

NU  • 

8 7 6 

5 4 

3 2 

1 

OEO 

3 

0.000 

0.000 

0.000 

0.000 

OEO 

3 

0.000 

0.000 

0.000 

0.008 

OEO 

3 

0.000 

0.016 

0.000 

0.008 

uEO 

= 

0.000 

0.016 

0.000 

0.008 

OEO 

a 

0.000 

0.016 

0.000 

0.008 

OEO 

3 

0.000 

0.016 

0.000 

0.008 

GEO 

» 

0.000 

0.016 

0.000 

0.008 

GEO 

a 

0.000 

0.016 

0.006 

0.043 

MA  » 

1.000E*-03  1.1 

300E+01 

HA  = 

1.250E-I 

01  1.1 

OOOE-03 

HA  =■ 

1.250E-I 

01  l.OOOE-03 

HA  ’ 

1.250E- 

01  1.1 

300E-03 

HA  " 

1.250E-I 

01  l.OOOE-03 

HA  ^ 

1 .250E-I 

01  1.1 

300E-03 

MA  =■ 

1.250E-' 

01  l.i 

300E-03 

HA  » 

4.480E*00  3.070E-02 

PCC 

3 

0.0000 

see 

= 0.0000 

PCC 

= 

0.0000 

see 

=■  0.0000 

PCC 

3 

0.0000 

see 

’ 0.0000. 

PCC 

3 

-0.3500 

see 

= 0.0000 

PCC 

3 

-0.3500 

see 

= 0.0000 

PCC 

= 

-0.3500 

see 

= 0.0000 

PCC 

= 

-0.3500 

see 

=■  0.0000 

PCC 

3 

-0.3500 

see 

= 0.0000 

PCC 

3 

-0.3500 

see 

= 0.0000 

PCC 

3 

0.0000 

see 

= 0.0000 

TO  = 0.00000  TS  » 0.00100  TE  * 0.30000 
MAXF  = 20  HAXM  =■  7 

NJTQS  =■  7 


MT 

1>  1.  1.  2. 

2.  0.750 

0.000  0.000 

0.000 

HT 

2»  1.  2.  3. 

1.  0.750 

0.000  0.000 

0.000 

MT 

3=  1.  3.  4. 

1.  0.750 

0.000  0.000 

0.000 

MT 

4=  1.  4.  5. 

1.  0.750 

0.000  0.000 

0.000 

MT 

5=  1.  5.  6. 

1.  0.750 

0.000  0.000 

0.000 

MT 

6»  1.  6.  7. 

1.  0.750 

0.000  0.000 

0.000 

MT 

7=  1.  7.  8. 

3.  0.750 

0.000  0.000 

0.000 

KFSM 

1 » 3 RFSM  = 

7 

FSM 

• 3.000E+00 

3.000E+00 

3.000E+00 

FSM 

» -l.OOOE+00 

-1.350E+00 

-6.500E-01 

FSM 

• 2.600E+02 

2.600E+02 

2.600E+02 

FSM 

• O.OOOE+00 

-3.500E-01 

3.500E-01 

FSM 

» O.OOOE+00 

O.OOOE+00 

0.000E*00 

FSM 

» l.OOOE+00 

6.500E-01 

1.350E+00 

FSM 

=■  -2.600E+02 

-2.600E+02 

-2.600E+02 

ACC 

= 1 0 

ACC 

» 0 2 

ACC 

* 0 2 

ACC 

« 0 2 

ACC 

=■  0 2 

ACC 

=>  0 2 

ACC 

=■  0 2 

ACC 

=■  0 2 

KFSA  =■  2 RFSA  = 

21 

FSA 

=■  l.OOOE+01 

4.000E<-00 

FSA 

» O.OOOE+00 

O.OOOE^OO 

FSA 

• O.OOOE+00 

O.OOOE+00 

FSA 

- 4.000E-02 

4.900E-02 

FSA 

» O.OOOE+00 

O.OOOE+00 

FSA 

• 8.200E-02 

5.000E-02 

FSA 

• -1.400E+02  ■ 

-9.810E+00 

FSA 

= l.OOOE-01 

1.000E<-00 

FSA 

» -l.lOOE+02  • 

-9.810E+00 

FSA 

» 1.060E-01 

O.OOOE+OO 

FSA 

= -1.900E+02 

O.OOOE+00 

FSA 

= 1.120E-01 

O.OOOE+00 

FSA 

» -1.900E+02 

O.OOOE+00 

FSA 

=■  1.150E-01 

0.OOOE*OO 

FSA 

» -1.250E+02 

O.OOOE+00 

FSA 

» 1.400E-01 

O.OOOE+00 

FSA 

» -1.500E+02 

O.OOOE+00 

FSA 

• 2.250E-01 

O.OOOE+00 

FSA 

» O.OOOE+00 

O.OOOE+00 

FSA 

• 2.800E-01 

O.OOOE+00 

FSA 

• O.OOOE+00 

O.OOUE^OO 

NET 

= 0 

N0T 

• 0 
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Reprint  of  input  data  of  second  CAD  model  (CAD2) 


RUN  2 


NST  = 

1 

ST  => 

8 1 

LBR  = 

3 

NU  =■ 

8 7 6 5 4 3 2 

1 

liEO  =■ 

0.000  0.000  0.000 

0.000 

GEO  =■ 

0.000  0.000  0.000 

0.008 

GEO  ’ 

0,000  0.016  0.000 

0.008 

GEO  • 

0.000  0.016  0.000 

0.008 

GEO  = 

0.000  0.016  0.000 

0.008 

GEO  = 

0.000  0.016  0.000 

0.008 

GEO  = 

0.000  0.016  0.000 

0.008 

GEO  = 

0.000  0.016  0.006 

0.043 

MA  > 

l.OOOE+03  l.OOOE+01 

MA  = 

1.250E-01  l.OOOE-03 

MA  > 

1.250E-01  l.OOOE-03 

MA  = 

1.250E-01  l.OOOE-03 

MA  => 

1.250E-01  l.OOOE-03 

MA  » 

1.250E-01  l.OOOE-03 

MA  ’ 

1.250E-01  l.OOOE-03 

MA  » 

4.480E+00  3.070E-02 

PCC 

" 0.0000 

see 

a 

0.0000 

PCC 

= 0.0000 

see 

= 

0.0000 

PCC 

= 0.0000 

see 

= 

0.0000 

PCC 

=■  -0.3500 

see 

a 

0.0000 

PCC 

= -0.3500 

see 

a 

0.0000 

PCC 

= -0.3500 

see 

a 

0.0000 

PCC 

= -0.3500 

see 

= 

0.0000 

PCC 

= -0.3500 

see 

3 

0.0000 

PCC 

=>  -0.3500 

see 

= 

0.0000 

PCC 

> 0.0000 

see 

= 

0.0000 

TO  =>  0,00000  TS  • O.OOlOO  TE  = 0.30000 
MAXF  =■  20  HAXM  = 7 

NJTQS  =■  7 


MT 

1=  1.  1.  2 

2.  0.750 

0.000 

0.000 

0,000 

MT 

2=  1.  2.  3 

1.  0.750 

0.000 

0.000 

0.000 

MT 

3=  1.  3.  4 

1.  0.750 

0.000 

0.000 

0,000 

MT 

4=  1,  4.  5 

1.  0.750 

0.000 

0.000 

0.000 

MT 

5=  1.  5.  6 

1.  0.750 

0.000 

0.000 

0.000 

MT 

6»  1.  6.  7 

1.  0.750 

0.000 

0.000 

0.000 

MT 

7=  1.  7.  8 

3.  0.750 

0.000 

0,000 

0.000 

KFSM  » 3 RFSM  » 

7 

FSM 

» 3.000E+00 

3.000E+00 

2,U00E+00 

FSM 

=>  -l.OOOE+00 

-1.350E+00 

O.OOOE+OO 

FSM 

=■  2.600E+02 

2.600E+02 

O.OOOE+OO 

FSM 

• O.OOOE+00 

-3.500E-01 

l.OOOE+OO 

FSM 

• O.OOOE+00 

O.OOOE+OO 

-l.OOOE+03 

FSM 

» l.OOOE+00 

6.500E-01 

O.OOOE+OO 

FSM 

» -2.600E+02 

-2.600E+02 

O.OOOE+OO 

ACC 

• 1 0 

ACC 

» 0 2 

ACC 

• 0 2 

ACC 

• 0 2 

ACC 

» 0 2 

ACC 

• 0 2 

ACC 

=■  0 2 

ACC 

•02 

KFSA  • 2 RFSA  =■ 

21 

FSA 

• l.OOOE+01 

4.000E+00 

FSA 

• O.OOOE+00 

O.OOOE+OO 

FSA 

• O.OOOE+00 

O.OOOE+OO 

FSA 

» 4.000E-02 

4.900E-02 

FSA 

• O.OOOE+00 

O.OOOE+OO 

FSA 

= 8.200E-02 

5.000E-02 

FSA 

• -1.400E'-02 

-9.810E+00 

FSA 

= l.OOOE-01 

l.OOOE+00 

FSA 

• -l.lOOE+02 

-9.810E+00 

FSA 

• 1.060E-01 

O.OOOE+OO 

FSA 

• -1.900E+02 

O.OOOE+OO 

FSA 

• 1.120E-01 

O.OOOE+OO 

FSA 

» -1.900E+02 

O.OOOE+OO 

FSA 

• 1.150E-01 

O.OOOE+OO 

FSA 

• -1.250E+02 

O.OOOE+OO 

FSA 

• 1.400E-01 

O.OOOE+OO 

FSA 

• -1.500E+02 

O.OOOE+OO 

FSA 

• 2.250E-01 

O.OOOE+OO 

FSA 

• O.OOOE+00 

O.OOOE+OO 

FSA 

• 2.800E-01 

O.OOOE+OO 

FSA 

• O.OOOE+00 

O.OOOE+OO 

NET 

• 0 

NBT 

• 0 

Reprint  of  input  data  of  third  CAD  model  (CAD3) 
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<UN  3 


NST  = 

1 

ST  =■ 

13  1 

L3R  =■ 

13 

NU  =■ 

13  12  11  10  9 

8 7 

6 5 

NU  > 

3 2 1 

liEO  » 

0.000  0.000 

0.000 

0.000 

GEO  = 

0.000  0.000 

0.000 

0.008 

GEO  • 

0.000  0.016 

0.000 

0.008 

GEO  = 

0.000  0.016 

0.000 

0.008 

GEO  = 

0.000  0.016 

0.000 

0.008 

GEO  = 

0.000  0.016 

0.000 

0.008 

GEO  = 

0.000  0.016 

0.000 

0.008 

GEO  • 

0.000  0.016 

0.000 

0.008 

GEO  » 

0.000  0.016 

0.000 

0.008 

GEO  =■ 

0.000  0.016 

0.000 

0.008 

GEO  • 

0.000  0.016 

0.000 

0.008 

GEO  = 

0.000  0.016 

0.000 

0.008 

GEO  =■ 

0.000  0.016 

-0.013 

-0.035 

MA  > 

l.OOOE+03  1. 

000E*01 

MA  » 

1.250E-01  1. 

OOOE-03 

MA  > 

1.250E-01  1. 

OOOE-03 

MA  > 

1.250E-01  1. 

OOOE-03 

MA  > 

1.250E-01  1. 

OOOE-03 

MA  > 

1.250E-01  1. 

OOOE-03 

MA  > 

1.250E-01  1. 

OOOE-03 

MA  • 

1.250E-01  1. 

OOOE-03 

MA  » 

1.250E-01  1. 

OOOE-03 

MA  > 

1.250E-01  1. 

OOOE-03 

MA  = 

1.250E-01  1. 

OOOE-03 

MA  = 

1.250E-01  1. 

OOOE-03 

MA  =■ 

4.480E+00  3. 

070E-02 

PCC 

= 0.0000 

see 

= 

0.0000 

PCC 

=■  0.0000 

see 

3 

0.0000 

PCC 

= 0.0000 

see 

= 

0.0000 

PCC 

= -0.3500 

see 

= 

0.0000 

PCC 

» -0.3500 

see 

3 

0.0000 

PCC 

= -0.3500 

see 

= 

0.0000 

PCC 

= -0.3500 

see 

3 

o.ouoo 

PCC 

= -0.3500 

see 

= 

0.0000 

PCC 

= -0.3500 

see 

= 

0.0000 

PCC 

=>  -0.3500 

see 

= 

0.0000 

PCC 

= -0.3500 

see 

= 

0.0000 

PCC 

=■  -0.3500 

see 

= 

0.0000 

PCC 

= -0.3500 

see 

3 

0.0000 

PCC 

» -0.3500 

see 

= 

0.0000 

PCC 

» 0.0000 

see 

= 

0.0000 

TO  ■ 

• 0.00000 

TS  = 

0. 

,00100  TE  =■  0.30000 

MAXF  » 20  MAXM  » 12 
NJTQS  » 12 


MT 

!=■ 

1. 

1. 

2. 

2. 

0.750 

0.000 

0.000 

0.000 

MT 

2» 

1. 

2. 

3. 

1. 

0.750 

0.000 

0.000 

0.000 

MT 

3» 

1. 

3. 

4. 

1. 

0.750 

0.000 

0.000 

0.000 

MT 

4» 

1. 

4. 

5. 

1. 

0.750 

0.000 

0.000 

0.000 

MT 

5= 

1. 

5. 

6. 

1. 

0.750 

0.000 

0.000 

0.000 

MT 

6* 

1. 

6. 

7. 

1. 

0.750 

0.000 

0.000 

0.000 

MT 

7> 

1. 

7. 

8. 

1. 

0.750 

0.000 

0.000 

0.000 

MT 

8» 

1. 

8. 

9. 

1. 

0.750 

0.000 

0.000 

0.000 

MT 

9= 

1. 

9. 

10. 

1. 

0.750 

0.000 

0.000 

0.000 

MT 

10» 

1.  10. 

11. 

1. 

0.750 

0.000 

0.000 

0.000 

MT 

11» 

1.  11. 

12. 

1. 

0.750 

0.000 

0.000 

0.000 

MT 

12» 

1.  12. 

13. 

3. 

0.750 

0.000 

0.000 

0.000 

KFSM  =» 

3 

RFSM 

= 

7 

FSM  = 

3 

FSM  =■ 

-1 

FSM  =■ 

2 

FSM  = 

0 

FSM  » 

0 

FSM  » 

1 

FSM  =■ 

-2 

ACC  =• 

1 

ACC  = 

0 

ACC  =« 

0 

ACC  = 

0 

ACC  = 

0 

ACC  =■ 

0 

ACC  = 

0 

ACC  =• 

0 

ACC  =■ 

0 

ACC  =■ 

0 

ACC  = 

0 

ACC  =■ 

0 

ACC  = 

0 

KFSA  » 

FSA  = 

1, 

FSA  =■ 

0, 

FSA  =• 

0. 

FSA  = 

4, 

FSA  = 

0, 

FSA  » 

8, 

FSA  = 

-1, 

FSA  =■ 

1. 

FSA  = 

-1, 

FSA  =■ 

1, 

FSA  =■ 

-1, 

FSA  =■ 

1, 

FSA  » 

-1 

FSA  » 

1 

FSA  « 

-1 

FSA  = 

1 

FSA  » 

-1 

FSA  » 

2 

FSA  = 

0 

FSA  = 

2 

FSA  =■ 

0 

NET  > 

0 

N8T  =■ 

0 

.000E*-00 

.OOOE+00 

.600E+02 

,000E+00 

.OOOE+00 

.OOOE+00 

.600E+02 

0 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

RFSA  = 
OOOE+01 
OOOE<-O0 
.OOOE+OO 
.OOOE-02 
.OOOE+00 
.200E-02 
.400E+02 
•OOOE-01 
.lOOE+02 
.060E-01 
.900E+02 
.120E-01 
.900E+02 
•150E-01 
.250E+02 
.400E-01 
.500E+02 
.250E-01 
,000E+00 
,800E-01 
.000E>00 


3.000E*-00 

-1.350E+00 

2,600E^-02 

-3.500E-01 

O.OOOE+00 

6.500E-01 

-2.600E<-02 


21 

4.000E+00 

O.OOOE+00 

O.OOOE^OO 

4.900E-02 

O.OOOE+00 

5.000E-02 

■9.810E+00 

l.OOOE+00 

9.810E+00 

O.OOOE+00 

O.OUOE+00 

O.OOOE+00 

O.OOOE+00 

O.OOOE+00 

O.OOOE+00 

O.OOOE+00 

O.OOOE+00 

O.OOOE+00 

O.OOOE+00 

O.OOOE+00 

O.OOOE+00 


3.000E>OJ 

-6.500E-J1 

2.600E+02 

3.500E-01 

O.OOOE+00 

1.350E+00 

-2.600E+02 


Reprint  of  input  data  of  fourth  CAD  model  (CAD4) 


i^UN  4 


NST  = 

1 

ST  = 

8 1 

LBR  =■ 

8 

NU  = 

8 7 1 

6 5 4 

3 2 

1 

GEO  » 

0.000 

0.000 

0.000 

0.000 

GEO  =■ 

0.000 

0.000 

0.000 

0.008 

GEO  = 

0.000 

0.016 

0.000 

0.008 

GEO  => 

0.000 

0.016 

0.000 

0.008 

GEO  =■ 

0.000 

0.016 

0.000 

0.008 
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A®ENDIX  C 

COMPARISON  OF  MASS  DISTRIBUTION  DATA  OF  THE  PART  572  DUMMY  AND  RECOMMENDATIONS 
FOR  A STANDARD  DATASET  FOR  CVS  MODELING* 


Jac  Wismans, 

Research  Institute  for  Road  Vehicles  TNO, 
Delft,  The  Netherlands 


Introducti on 


Crash  Victim  Simulations  (CVS)  models  most  frequently  used  now  in  the  field  of 
automotive  safety  are  MVMA^^\  CALSPAN^^^  and  MADYMO^^^.  Datasets  used  in 
these  computer  programs  to  represent  the  50%  Part  572  dummy  appear  to  be  not 
completely  identical.  These  differences  are  mainly  due  to  differences  in: 

- segment  division 

- number  of  segments  used  to  describe  the  dummy 

- location  of  segment  coordinate  systems 

- methods  used  to  measure  dummy  characteri sti cs 

- the  dummies  selected  for  the  measurement  of  the 
dummy  characteri sti cs 

- program  features  like  e.g.  for  the  description  of  joint 
characteri sties 

In  addition  some  dummy  characteristics  are  hard  to  measure  or  cannot  be 
measured  at  all.  Also  the  type  of  simulation,  i.e.  frontal  versus  side  or 
pedestrian  impact,  restrained  versus  unrestrained,  2D  versus  3D  etc.  can 
effect  the  representation  of  the  dummy. 


* This  document  has  been  prepared  for  the  Analytical  Human  Simulation  Task 
Force. of  the  SAE  Human  Biomechanics  and  Simulation  Subcommittee  (aug.  1983). 
Comments  and  suggestions  can  be  sent  to  the  chairman  of  this  Task  Force: 
dr.  Priyaranjan  Prasad,  Ford  Motor  Company,  Automotive  Safety  Center, 

P.O.  Box  253,  Dearborn,  Mi  48221,  USA. 
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An  analysis  will  be  presented  here  of  Part  572  dummy  mass  distribution  data 
based  on  measurements  conducted  at  three  different  laboratories: 

- "Hubbard"  presented  in  Reference  (4). 

- "Calspan"  presented  in  Volume  2 of  Reference  (2). 

- "TNO"  presented  in  Reference  (5). 

"Hubbard"  data  don't  include  3-dimensional  moments  of  inertia  and  are  used  in 

datasets  for  frontal  impacts  of  MVMA  and  the  MADYMO  2-dimensional  version. 
"Calspan"  data  are  used  in  the  CALSPAN  model  and  "TNO"  data  in  the  MADYMO 
3-dimensional  version.  Mass  distribution  data  to  be  presented  here  include 
masses,  moments  of  inertia,  center  of  gravity  location  and  distances  between 
pi vots . 

On  the  basis  of  this  analysis,  a standard  mass  distribution  dataset  will  be 
proposed  derived  from  these  three  sources.  This  dataset  will  be  extended  in 
future  with  joint  characteristics  and  body  stiffness  data. 

Division  in  Segments 

In  order  to  measure  mass  distribution  data,  the  dummy  has  to  be  divided  into  a 
number  of  segments.  The  segment  division  used  by  the  three  laboratories  is  not 
completely  identical  which  makes  it  more  difficult  to  compare  in  a direct  way 
the  results  of  these  laboratories.  In  this  study,  a division  in  13  segments  is 
proposed:  head,  neck,  upper  torso,  center  torso,  lower  torso,  upper  arms, 
lower  arms  and  hands,  upper  legs,  lower  legs  and  feet.  The  center  torso  has 
2 parts:lumbar  spine  and  abdominal  insert.  Figure  1 shows  roughly  which 
components  belong  to  a segment.  Bolts  and  other  small  components  are  not 
separately  assigned  to  a segment  here,  since  difference  methods  are  used  in 
the  three  studies.  The  effects  of  this  is  expected  to  be  small.  Further  it  is 
assumed  that  no  instrumentation  is  included  in  the  dummy  segments.  Table  1 
summarizes  the  deviations  between  the  segment  division  used  in  this  study  and 
the  ones  used  by  the  laboratories. 

The  standard  dataset  for  a 3-dimensiona1  representation  of  the  dummy  in  a 
mathematical  model,  proposed  here,  will  be  based  on  a similar  division  in  13 


segments.  A corresponding  dataset  for  a 2-dimensional  representation  simply 
can  be  derived  from  this  proposed  3-dimensional  dataset  and  would  consist  of 
13  segments  in  case  of  side  impacts,  and  9 segments  in  case  of  frontal 
impacts . 

Segment  coordinate  systems 


For  each  segment,  a right-handed  coordinate  system  x,  y,  z is  defined.  In  a 
standing  upright  position  with  the  arms  down  and  the  legs  straight,  all  z-axis 
are  in  the  upward  direction,  all  x-axis  are  forward  and  all  y-axis  are  to  the 
left  (see  Figure  2).  The  location  of  the  local  coordinate  systems  is  indicated 
in  Figure  1 and  specified  in  more  detail  in  Table  2. 

The  origins  of  the  coordinate  systems  are  selected  in  joint  centers  of 
rotation  as  far  as  the  limbs  are  concerned  and  in  the  center  of  the  endplanes 
of  the  rubber  cylinders  representing  neck  and  spine,  for  the  other  segments. 

Since  in  the  dummy  the  neck  centerline  is  rotated  20  degrees  forward  relative 
to  the  spine  centerline  (2),  the  head  and  neck  local  coordinate  systems  are 
rotated  20  degrees  with  respect  to  the  upper  torso  cordinate  system  (if  no 
forces  are  applied  on  the  head/neck  structure).  In  the  reference  position 
(standing)  illustrated  in  Figure  2,  all  other  segment  coordinate  systems  have 
the  same  orientation. 

Mass  distribution  data 


In  this  section,  data  from  the  three  sources  are  compared  with  each  other. 
Where  possible,  additional  calculations  are  made  in  order  to  have  the  data 
expressed  in  agreement  with  the  preceding  segment  and  coordinate  system 
def i ni ti ons . 

Both  CALSPAN  and  TNO  data  are  based  on  measurements  of  one  dummy  (Alderson) 
while  the  data  presented  by  Hubbard  relate  to  two  dummies  by  different 
manufacturers  (Alderson  and  Sierra).  Average  values  will  be  presented  here  if 
separate  measurements  for  the  left  and  right  segments  were  conducted. 
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Table  3 shows  a comparison  of  measured  segment  masses.  The  largest  difference 
appears  to  be  in  the  lower  torso  mass  which  is  1.6  kg  higher  in  the  Sierra 
dummy  (Hubbard)  than  in  the  Alderson  dummy  (TNO).  Differences  in  mass  for  the 
other  segments  are  smaller  than  0.5  kg. 

Results  for  the  center  of  gravity  measurements  are  shown  in  Table  4.  Values 
between  -5  and  +5  mm  are  neglected  in  this  table.  It  turns  out  that  except  for 
the  upper  leg,  all  centers  of  gravity  are  located  in  the  (x,z)-plane  of  the 
segment.  It  should  be  noted  that  results  of  the  CALSPAN  measurements  for  the 
lower  torso  and  the  upper  leg  cannot  be  compared  with  the  other  sources  due  to 
different  segment  definitions.  Differences  in  center  . of  gravity  locations 
between  various  measurements  are  found  to  be  relatively  small.  Largest 
differences  (i.e.  3 cm)  appear  to  be  in  the  x-coordinate  of  the  abdominal 
i nsert . 

Table  5 summarizes  the  segment  principal  moments  of  inertia.  Only  Hubbard 
conducted  measurements  about  the  local  y-axis.  In  the  TNO  measurements,  the 
axis  of  the  local  coordinate  systems  were  assumed  to  have  the  same  orientation 
as  the  principal  axis.  According  to  the  measurements  of  CALSPAN,  this  appears 
to  be  a reasonable  approximation  for  most  of  the  dummy  segments.  Significant 
differences  between  both  axis  systems  were  noted  only  for  the  head  and  the 
lower  torso  as  indicated  in  Table  5 (note  4 and  note  7). 

Like  for  the  center  of  gravity  measurements,  CALSPAN  moments  of  inertia  for 
the  lower  torso  and  upper  leg  cannot  be  compared  with  the  results  of  the 
other  laboratories  due  to  differences  in  segment  definition.  The  same  holds, 
for  the  upper  torso  measurements  by  TNO,  which  were  based  on  another  segment 
definition  (i.e.  without  shoulders).  The  deviations  between  moment  of  inertia 
are,  in  general,  less  than  20%,  except  for  smaller  segments  (i.e.  spine  and 
neck)  and  for  moments  of  inertia  about  the  longitudinal  axis  of  slender 
elements  (limbs).  In  the  case  of  the  neck  and  spine,  calculated  values 
(CALSPAN/TNO)  are  expected  to  be  more  accurate  than  measured  ones  (Hubbard). 

Table  6 presents  the  positions  of  the  origins  of  the  segment  coordinate 
systems  expressed  in  the  coordinates  of  an  adjoining  element.  For  the  limbs, 
the  origins  are  taken  identical  to  joint  centers  of  rotation.  For  the  torso. 
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head  and  neck  segments,  no  unique  pivots  can  be  defined  due  to  the  flexible 
spine  and  neck  structures.  It  will  be  ass^Jmed  here  that  these  joints  can  be 
reasonably  located  on  the  centerlines  in  the  endplanes  of  the  rubber  cylinders 
so  coinciding  with  origins  of  torso,  neck  and  head  segments.  Thus,  this  table 
implicitly  presents  the  joint  locations  (and  link  length).  Data  provided  by 
TNO  are  based  on  measurements  from  the  assembly  drawings,  so  are  expected  to 
be  less  accurate.  Measured  data  by  Calspan  and  Hubbard  only  show  significant 
differences  for  the  shoulder  (0.03  m)  and  for  the  knee  (0.01  m).  With  respect 
to  the  shoulder,  this  most  likely  will  be  due  to  differences  in  defining  the 
reference  position  of  the  shoulder  with  respect  to  the  torso. 

Recommendation  for  a dataset  (mass  distribution) 


Based  on  the  preceding  analysis,  a dataset  will  be  proposed  here  which  can  be 
used  for  3-dimensional  CVS  models.  This  dataset  is  presented  in  Table  7 (mass 
and  moments  of  inertia).  Table  8 (center  of  gravity)  and  Table  9 (joint 
location).  The  dataset  only  holds  for  a division  of  the  dummy  in  13  segments 
and  for  a segment  and  coordinate  system  definition  as  proposed  in  this  study 
(see  Figure  1-2  and  Table  1,  2).  Footnotes  in  the  tables  describe  in  what  way 
the  data  were  obtained. 

In  this  way  the  complete  mass  distribution  of  the  Part  572  dummy  is  defined. 
It  should  be  noted  here  that: 

- in  a standing  position  of  the  dummy  (see  Figure  2)  all  local 
coordinates  except  for  the  head  and  the  neck  have  the  same  orientation. 
The  neck  and  head  local  coordinates  system  are  rotated  forward  20°  with 
respect  to  the  torso  local  coordinate  system. 

- all  local  principal  axis  except  for  the  head  are  assumed  to  have  the 
same  orientation  as  the  proposed  local  coordinate  systems.  For  the  head 
the  principal  axis  system  is  taken  -43  degrees  rotated  about  the  y-axis 
of  the  local  head  coordinate  system  (i.e.  rotated  backwards). 
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Table  1 Division  of  Part  572  Dummy  in  13  Segments 


NAME 

lower  torso 


center  torso 


upper  torso 

neck 

head 

left  (right)  upper  arm 

left  (right)  lower  arm 
and  hand 

left  (right)  upper  leg 

left  (right)  lower  leg 
and  foot 


COMMENTS 

CALSPAN;  abdominal  insert,  lumbar 
spine  and  femur  and  retainer 
assembly  E05/E06  are  incorporated 
in  the  lower  torso  segment 

Consists  of  2 separate  segments: 
lumbar  spine  and  abdominal  insert. 
Data  for  both  components  will  be 
given  separately. 

Includes  shoulder.  TNO  and 
Hubbard:  separate  measurements  for 
the  shoulder. 


Hubbard:  separate  measurements  for 
the  hand  and  lower  arm. 

CALSPAN:  without  femur  and 

retainer  assembly  E05/E06. 

Hubbard  and  CALSPAN:  separate 
measurements  for  lower  leg  and 
foot . 
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Table  2 Specification  of  local  coordinate  system*  (see  Figure  1) 


NAME 

COMMENTS 

lower  torso/central  torso 

origin: 

centerline  of  lumbar  spine  at 
spine/lower  torso  attachment 
surface. 

z-axis : 

along  centerline  lumbar  spine 

upper  torso 

origi n : 

centerline  of  lumbar  spine  at 
spine/upper  torso  attachment 
surface  (i.e.  in  plane  of 
bottom  side  of  D/03  lumbar 
adapter ) . 

z-axi s : 

along  centerline  lumbar  spine. 

neck 

ori  gi  n : 

centerline  of  neck  at  upper 
torso/neck  attachment  surface. 

z-axi s : 

along  centerline  neck. 

head 

origin: 

centerline  of  neck  at  head/neck 
attachment  surface. 

z-axi s : 

along  centerline  neck. 

upper  arm 

ori gi n : 

center  of  shoulder  pivot  for 
abduction/adduction  motions. 

x-axi s : 

along  centerline  of  shoulder 
pivot  for  abduction/adduction 
moti ons . 

z-axis: 

along  the  line  between  center  of 
shoulder  and  elbow  pivot. 

lower  arm 

origi  n : 

center  of  elbow  pivot. 

y-axis : 

along  centerline  of  elbow 

z-axi s : 

along  the  line  between  the  pivot 
center  of  the  elbow  and  the 
center  of  the  wrist  pivot. 

upper  leg 

ori  gin: 

center  of  the  hip. 

y-axi  s : 

parallel  to  centerline  of 

knee-pivot  bolt. 

z-axi s : 

intersects  the  knee  centerline. 

lower  leg 

origin: 

center  of  knee. 

y-axi s : 

along  the  centerline  of  the  knee 

. 

pivot  bolt. 

z-axi s : 

along  a line  from  the  center  of 
the  ankle  joint  (flexion  joint) 
to  the  center  of  the  knee  joint. 

* Note:  if  not  specified  in  another  way,  all  positive  z-axis  are 
upward,  all  positive  x-axis  are  forward  and  all  positive  y-axis 
are  to  the  left  (in  standing  position). 


Table  3 Segment  Masses 


NAME 

SIERRA 

HUBBARD 

(kg) 

HUBBARD 

(kg) 

ALDERSON 

CALSPAN 

(kg) 

TNO 

(kg) 

lower  torso 

12.7 

11.6 

11.62^^ 

11.1 

lumbar  spine 

1.4 

1.29 

1.36 

1.33 

abdominal  insert 

1.2 

1.4 

1.44 

1.32 

2) 

upper  torso 

17.4 

17.7 

17.16 

17.18 

neck 

0.9 

0.9 

0.82 

0.89 

head 

4.5 

4.6 

4.38 

4.2 

one  upper  arm 

2.2 

2.2 

2.16 

2.3 

one  lower  arm  & one 

hand 

2.15 

2.2 

2.09 

2.16 

one  upper  leg 

9.9 

9.85 

9.55^> 

9.41 

one  lower  leg  & one 

foot 

4.4 

4.3 

4.42 

4.54 

4) 

total  dummy  mass 

75.4 

74.59 

73.22 

72.84 

1)  After  correction  for  mass  abdominal  insert  (-  1.44  kg),  lumbar 
spine  (-  1.36  kg),  and  femur  and  retainer  assembly  (-  1.4  kg). 

2)  Includes  shoulders  (mass  of  one  shoulder  is  1.6  kg  (Alderson)  and 
1.5  kg  (Sierra)  according  to  Hubbard  measurements  and  1.45  kg 
according  to  TNO  measurements). 

3)  After  correction  for  mass  of  femur  and  retainer  assembly  (+  1.4  kg) 

4)  Sum  of  separate  segments. 


lable  ^ 

Center  of  Gravity 
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TfJO  arid  HoLitard.  average  of  right  and  left  segments;  CALSPAN:  measurement  of  rinht  seoment.  Coordinates  presented  in  the  coordinate 

system  of  a left  segment 
fJot  measured. 


Table  7 Model  Dataset:  Masses  and  Moments  of  Inertia 


1) 


NAME 

MASS^^ 

Moments  of  Inertia 

kg 

V 

Ix  kgm 

ly  kgm'^ 

Iz  kgm 

1 ower 

torso 

11.76 

0.1297^^ 

0.0817^^ 

0.1393^^ 

center  torso 

2.69 

0.0140^^ 

0.0159^^ 

0.0186^^ 

upper 

torso 

17.36 

0.2352^^ 

0.1896^^ 

0.1508^^ 

neck 

0.88 

0.0014^^ 

0.0014^^ 

0.0006^^ 

head 

4.42 

0.0248^^ 

0.0307^^ 

0.0184^^ 

upper 

arm 

2.22 

0.0161^^ 

0.0156^^ 

0.0016^^ 

1 ower 

arm 

2.15 

0.0311^^ 

0.0301^^ 

0.0014^^ 

upper 

leg 

9.68 

0.1300^^ 

0.1387^^ 

0.0170^^ 

1 ower 

leg 

4.42 

0.1315^^ 

0.1271^^ 

0.0046^^ 

total 

dummy  mass 

74.05 

1)  Except  for  the  head,  orientation  of  principal  axes  is  taken 
identical  to  local  coordinate  axes.  For  the  head,  the  principal 
axis  system  is  taken  -43  degrees  rotated  about  the  y-axis  of  the 
local  head  coordinate  system  (i.e.  rotated  backward). 

2)  Average  Hubbard/CALSPAN/TNO 

3)  TNO 

4)  CALSPAN 

5)  Average  CALSPAN/TNO 

6)  Average  Hubbard/TNO 

7)  Average  Hubbard/CALSPAN. 


Table  8 Model  Dataset:  Location  Center  of  Gravity 


NAjME 

X (M) 

y (M) 

z (M) 

lower  torso^^ 

0.026 

0 

-0.079 

3) 

center  torso  ' 

0.033 

0 

0.072 

2) 

upper  torso  ' 

0.029 

0 

0.162 

neck^^ 

0 

0 

0.063 

head'^^ 

0.006 

0 

0.028 

2) 

upper  arm  ' 

0 

0 

-0.122 

2 ) 

lower  arm  ' 

0 

0 

-0.167 

(left)  upper  leg^^ 

0 

0.006 

-0.207 

2 ) 

lower  leg 

0.016 

0 

-0.272 

1)  Average  of  Hubbard/TNO 

2)  Average  of  Hubbard/CALSPAN/TNO 

3)  Calculated  from  lumbar  spine  and  abdominal 
insert  values  of  Hubbard/TNO 

4)  Average  of  Hubbard/CALSPAN. 


Table  9 Model  Dataset:  Location  Origin  of  Joints 


NAME 

Relative  to  coordinate 
system  of: 

X 

(m) 

y 

(m) 

z 

(m) 

lower  spine 

lower  torso 

0 

0 

0 

upper  spine^^ 

center  torso 

0 

0 

0.132 

lower  neck^^ 

upper  torso 

0.065 

0 

0.318 

upper  neck^^ 

neck 

0 

0 

0.124 

(left)  shoulder^^ 

upper  torso 

0.030 

0.189 

0.260 

el bow^  ^ 

upper  arm 

0 

0 

-0.261 

(left)  hip^^ 

lower  torso 

0.042 

0.087 

-0.072 

(left)  knee^^ 

upper  leg 

0 

0.008 

-0.405 

Ij  Average  of  Hubbard/Cal span 
2)  Calspan. 
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